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ABSTRACT 


The effects of power system faults upon a generator with a super- 
conducting field winding are presented and discussed. Induced current, 
heating, and electromechanical stress effects upon the superconducting 
field winding and rotor structure are investigated. Design criteria 
are developed which insure that a superconducting generator can survive 
a power system fault. 


Two classes of faults are considered and design criteria relating 
to each are developed. Survival of a short-term fault which is cleared 
by a system circuit breaker requires maintenance of super conductivity; 
survival of a sustained fault which cannot be cleared requires mainte- 
nance of machine structural integrity without permanent damage. 


An example design for a 1000 MVA superconducting generator is 
developed. The design illustrates that it is possible to satisfy the 
criteria for fault survival, while maintaining the advantages of lower 
weight, lower synchronous reactance, improved efficiency, better damp- 
ing characteristics, larger iets and higher rotor natural frequencies 
than for a comparable conventional generator. 
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Glossary of Symbols (continued) 
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CHAPTER I 


Introduction and Statement of the Problem 


The increasing demand for electric power has made practical the use 
of higher and higher rated generating units. Economic advantages are 
realized in the capital and operational costs of both the prime mover 
and the generator as the unit rating increases [1]. A point of diminishing 
returns in economic advantage is reached for the generator above ratings 
of about 1000 MVA for the reasons given below. Generators with super- 
conducting field windings can extend this limit; thus, the economies of 
scale can be realized in even higher rated units. [2] 

The economy realized in conventional generators from increased 
rating is the result of better armature cooling and thus more effective 
use of the armature copper. However, as the rated current for a given 
machine is increased, the per-unit reactances of the machine are increased 
and the inertia constant is decreased, both of which degrade the transient 
stability margin. To decrease synchronous reactance, air gap lengths 
must be increased, resulting in the requirement for more excitation power 
and greater power loss in the field winding. Hence, the larger rated 
machine with the better utilized armature becomes less efficient and 
less economical in an operational sense. The crossover occurs at a 
rating of about 1000 MVA. 

The capital cost advantage associated with higher rated units is 
offset above 1000 MVA for conventional generators by construction and 
shipping considerations. The stator iron makes the larger machines so 


heavy that they must be shipped in pieces and assembled at the site. 
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Rotor diameters are limited by mechanical stresses requiring length 
increases for increased power ratings. Therefore, weight increases pro- 
portional to length further complicate shipping and assembly. Lower rotor 
mechanical natural frequencies are another result of increased length. 

Generators with superconducting field windings do not require iron 
in the rotor or stator to produce high magnetic flux density at the 
armature bars. With no resistance loss in the field winding, and no 
concern with iron saturation, the ampere turns of the field winding can be 
increased to produce up to about 50 kilogauss in the machine center and 
20 to 30 kilogauss in the armature region. 

Elimination of iron has several beneficial effects. First, machine 
weight is reduced; also, the elimination of iron within the armature pro- 
vides space for more armature conductors. The machine volume for a 
given power rating is thus reduced. Insulation problems in the armature 
are also reduced because of the absence of iron at ground potential; 
machines can thus be designed for higher voltages. Second, machines 
without iron have lower reactances. Even though the inertia is reduced, 
transient stability margin has been shown to improve over conventional 
machines [5]. Third, the superconducting generator is slightly more 
efficient because of the elimination of field winding resistance losses. 

The configuration for a superconducting generator is as shown in 
Fig. I-l. The armature is at room temperature; it is similar to a con- 
ventional armature except for the absence of iron interleaved with the 
conductor bars. A winding scheme developed by Kirtley and Smith [21] 
permits use of much less insulation, and provides better utilization of 


the armature space for carrying current. A laminated iron shield 
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surrounds the armature to provide a uniform boundary for the magnetic 
fields and eliminate unbalanced forces upon the field winding. 

The superconducting field winding is located inside a liquid helium 
space in the rotor and is maintained at 4.2° K. Surrounding the field 
winding is a thermal shield maintained at approximately 20° K to inter- 
cept thermal radiation from the room temperature parts. The outer part 
of the rotor is a conducting electrical shield to intercept asynchronous 
magnetic fields produced in the armature. If the superconducting field 
winding were exposed to these asynchronous fluxes, losses would be pro- 
duced and the field winding could be driven normal. The layer of stainless 
steel beneath the electrical shield provides structural support for the 
shield when it is subjected to large electromechanical stresses under 
fault conditions. 

Superconducting generators have been shown to possess several 
operational and economic advantages over conventional machines. However, 
the effects of faults upon these generators have not been studied. The 
large oscillating torques and large armature currents producing rotor 
heating following a short circuit on a power system have been known and 
studied for many years. The internal effects produced by these fault 
conditions would seem to be more severe in a superconducting generator 
for the following reasons: 

1) The thermal isolation of the rotor to minimize conduction heat 
leak into the cryogenic region makes it structurally more 
vulnerable to large electromechanical stresses. 

2) Low thermal capacity of rotor parts at cryogenic temperatures 


can result in a temperature rise due to rotor heating during 
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faults. Small temperature rise in the superconductor can cause 
it to go normal. 
3) Induced currents in the field winding may cause it to go normal 

if the critical current is reached, or if hysteresis losses raise 

its temperature sufficiently. 
The first consideration in making superconducting generators able to 
withstand the effects of faults is the superconducting winding itself. 
If a short-term fault drives the superconductor normal, it might be neces- 
Sary to remove the generator from the system, whereas maintenance of 
superconductivity would probably make continued operation possible after 
the fault is cleared. The superconductor can be driven normal if (a): 
the current rises above the critical value; (b): the temperature rises 
due to heat transfer from other rotor parts, or due to internal heat 
generation in the superconductor itself; (c): possibly, excess motion of 
the superconductor results during the fault. 

The electrical shields, cylindrical conducting shells around the 
field winding, serve to mitigate all of these fault effects. A large 
portion of the induced currents, eddy current heating, and oscillating 
stresses associated with the period following a fault, are transferred 
from the superconducting winding to these shields. The problems of 
keeping the superconductor superconducting are, in part, traded off for 
problems of making the electrical shields able to withstand the fault 
effects. 

The second consideration in making the superconducting generator 


able to survive faults concerns the structural integrity of the rotor. 


a a 





17 


The generator rotor is subjected to large electromechanical stresses 
following a fault. These stresses include not only the well known shear 
Stresses associated with fault torques, but also normal stresses, because 
of the initial exclusion of demagnetizing armature flux from the rotor. 
The rotor torque~supporting structure at the rotor ends of a superconduct- 
ing generator is minimized to reduce heat leak by conduction to the cryo- 
genic region. In addition, much of the fault torque is taken by the 
electrical shield, which is of relatively thin material. Hence, both 

the rotor thermal distance pieces and the electrical shield are vulner- 
able to fault torques. The electrical shield is also subjected to the 
normal electromechanical stresses. Its only radial support is at the 
ends, to prevent thermal communication directly to the field winding 
repion, making it particularly susceptible to deformation from these 
radial stresses. 

The electrical shield has evolved from the standard thermal radia- 
tion shield used in helium dewars. The original function of the thermal 
shield was to intercept steady~state radiant energy from the room-temp- 
erature stator. Also, it intercepts electromagnetic energy in the form 
of pe enaeaots mapnetic fields produced by armature phase imbalance, 
and serves as a damper winding to damp machine oscillations. This study 
shows that it can also provide protection for the superconducting wind- 
ing during fault conditions.Providing protection from asynchronous mag- 
netic fields is shown to imply that the shield intercepts the transient 


electromechanical stresses during faults. 
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To provide sufficient mechanical support for the electrical shield, 
it is more efficient to place the electrical shield at room temperature 
to minimize conduction heat leak into the crvogenic region. The function 
of thermal radiation shielding is then transferred to a radiation shield 
at cryogenic temperatures. The electrothermal shield proposed by Thullen 
[17] has thus become an electrical shield at room temperature to inter- 
cept asynchronous flux and take the electromechanical stresses, and a 
thermal shield at cryogenic temperature to intercept thermal radiation. 

This study of fault effects in the superconducting generator reveals 
additional benefits of the electrical and thermal shields for fault pro- 
tection, but at the same time yields requirements upon their desipn to 
insure the ability of the generator to survive faults. This study is an 
attempt to develop the analytical tools necessary to predict the environ- 
ment to which the rotor is subjected following the fault. The electrical, 
thermal, and mechanical responses of the affected parts are investigated 
in an attempt to identify potential critical failures. The steady-state 
aspects of the shields’ performances are also compiled. Sample calcula- 
tions have been made for the existing 80 KVA machine, the 2 MVA machine 
now being built, and power system size generators of lLOOOMVA rating, 
which is included in Chapter II to indicate the relative importance of 
the fault effects and the measures necessarv to protect against them. 

The ultimate test of how a superconducting generator performs under 
fault conditions can be answered only with some experimental experience 
with a real generator subjected to a fault. This study is primarily 
analytical, but supporting experiments to check critical aspects of the 


analytical models were performed and are reported in Chapter V. 
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The existing 80 KVA generator was subjected to a three-phase fault at 
low level. More experiments are planned for this machine and for the 


2 MVA macnine being built. 
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CHAPTER II 
Results and Example Design 

Tne purpose of this study is to evaluate the performance of a power- 
system size superconducting generator under fault conditions normally 
encountered in operation. The design requirements necessary to insure 
that the generator can survive fault conditions are to be identified and 
established. Two classes of fault conditions will be considered. The 
first will be a short-term fault which is cleared by a system circuit 
breaker. The duration of a short-term fault will be five to 15 cycles, 
depending on the type of circuit breaker. The second class of fault will 
be a sustained fault which occurs on the machine side of the system cir- 
cuit breaker, such as the machine terminals, and cannot be cleared by 
circuit breaker action. 

The criteria for survival of these two classes of faults will be as 
follows. To survive the short-term fault, superconductivity of the field 
winding must be maintained. If the superconducting winding is driven 
normal during the period before the fault is cleared, the machine must be 
removed from the system. If superconductivity is maintained, the machine 
can continue to supply power to the system through the remaining unfaulted 
lines after the fault has been cleared. To survive a sustained fault, 
structural integrity of the generator must be maintained so that perma- 
nent damage does not result. 

Maintenance of superconductivity through a short-term fault requires 
that heat dissipation be limited in the cryogenic region where the field 
winding is located, and that the level of current induced in the field 


winding be limited. Superconductors can be driven normal, if their 
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temperature rises more than a few degrees, or if their current level 
increases above a critical level. Maintenance of structural integrity 
during a sustained fault requires that the rotor parts subjected to the 
large electromechanical stresses be substantial enough not to yield. The 
yielding criterion also applies to the short-term fault, but it is more 
critical for the sustained fault because of the large amount of heating 
due to the dissipation from induced currents in the rotor structure. 

To maintain superconductivity, a rotating, electrically conducting 
shield must be provided to intercept the asynchronous magnetic fields pro- 
duced by trapped armature fluxes and unbalanced conditions of the arma- 
ture phases. The damper windings or bars of conventional generators serve 
this purpose, and limit heating of the rotor iron. However, for the super- 
conducting generator, the electrical shield must not have thermal comnu- 
nication with the field winding region. The induced currents and dissi- 
pation must be intercepted outside the cryogenic region. This study shows 
that a shield which intercepts the induced currents and dissipation also 
intercepts the large electromechanical stresses associated with faults. 

To support these stresses, considerable mechanical structure must be pro- 
vided for the electrical shield. A room-temperature electrical shield 

appears to be the most economical design, because of the large conduction 
heat leak through the structure which would result for a cryogenic shield. 

To intercept thermal radiation from the room-temperature electrical 
shield, a secondary thermal shield is provided at a temperature somewhat 
above that of the field winding (e.g., 20° K). The electrothermal shield 
proposed by Thullen [17] thus becomes a room-temperature electrical shield 


to intercept dissipation and electromechanical stresses and a cryogenic 
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thermal shield to intercept thermal radiation. 

The transfer of transient electromechanical stresses from field wind- 
ing to a room-temperature electrical shield is very desirable, since it 
permits the use of minimum structure to support the field winding, and 
thus limit heat conduction into the 4.2° K region. However, during fault 
conditions, the power dissipation rate in the electrical shield far 
exceeds the cooling rate; therefore the electrical shield's temperature 
will rise, depending on its thermal capacity. Since structural proper- 
ties of metals are degraded as their temperatures rise above room tempera- 
ture, the ability of the electrical shield to withstand the electromechan- 
ical stresses is related to its thermal capacity. 

In addition to the fault survival criteria which affect the electri- 
cal shield design, certain steady-state requirements must be considered. 
First, the electrical shield will serve as a damper winding to damp 
machine oscillations following a transient. Second, a steady-state elec- 
trical dissipation will occur due to slight armature current phase unbal- 
ance. The rotor cooling system must be designed to remove this heat. To 
minimize steady-state losses in the electrical shield, a highly electric- 
ally conducting material is desirable. This also makes the shield effec- 
tive in intercepting asynchronous magnetic fields during fault conditions. 
But a highly conducting shield with sufficient thickness to provide ther- 
mal capacity provides very little damping. Conflicting requirements 
therefore result from shielding,steady-state dissipating, and damping. 

The requirement upon the electrical shield design can be summarized 


as follows: 
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1) Keep steady-state power dissipation to a level acceptable for 
a reasonable cooling systen. 

2) Provide as much damping of machine oscillations as possible. 

3) Maintain sufficient shielding of asynchronous magnetic fields 
to protect the superconducting winding in the steady state, and 
during short-term faults. 

4) Provide thermal capacity to absorb transient dissipation without 
excessive temperature rise for a sustained fault. 

5) Provide sufficient structure to support the transient electro- 
mechanical stresses. 

6) Insure that no rotor mechanical natural frequencies are near 


frequencies of electromechanical or rotational stresses. 


This study provides a quantitative expression of the above effects in 


order to arrive at an acceptable design which compromises between the 
conflicting requirements. The following sections of this chapter discuss 
the tradeoffs of shielding, damping, structural integrity, thermal capa- 
city, and steady-state performance. Section II-4 presents a design for a 
1000 MVA machine which is based upon these requirements. Chapters III 
and IV present in more detail the analytical tools which have been devel- 
oped to evaluate the fault effects and steady-state requirements on the 
superconducting generator. Chapter V is a summary of experiments which 
have been performed to lend credence to some of the analytical expres- 
Sions. Chapter VI is a summary of results and conclusions with suggest- 


ions for further work. 
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II-1 Keeping the Superconductor Superconducting 


During the period immediately following a fault, the most adverse 
conditions are developed, which can tend to drive the superconducting 
field winding normal. Alternating currents are induced in the field 
winding, as well as a rise in the de current level due to the demagnetizing 
effect of the armature fault currents. A large heat load is imposed in 
the rotor, due to trapped flux in the armature, and resulting imbalance of 
shorted armature phases. Alternating stresses are exerted in the rotor 
which can produce motion of the winding, and possibly degrade its per- 
formance. Just how critical each of these effects is in producing a super- 
conducting to normal conductivity transition in the winding will require 
experimental investigation, as can be done with the generator presently 
being built. This section is intended to illustrate the design require- 
ments for a shield to mitigate the adverse effects upon the field winding. 

The degree of protection which the electrical shield affords the 
superconducting field winding is related to the shield attenuation factor. 
This factor indicates the degree to which asynchronous magnetic flux, pro- 
duced in the stator, is excluded from the rotor, and is numerically equal 
to the ratio of the asynchronous magnetic field inside the shield to the 
asynchronous field which would exist without the shield. The factor is 
involved in expressions for induced field current, rotor heating in the 
cryogenic region, and mechanical stresses excited upon the field winding. 

The attenuation produced by an electrical shield is a steady-state 
phenomenon dependent on the frequency of the asynchronous magnetic flux. 
Associated with any conducting cylindrical shield there is an L/R time 


constant equivalent to the time constant of the damper winding circuit 
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used to model generators. Due to the process of magnetic diffusion in 
the continuum conducting shield, there are in fact an infinite number of 
such time constants, but the longest of these we will call the "electrical 


shield time constant", TS, and it is given by 


To - CAOUR 
A #= shield thickness 
R = shield radius 
o = shield electrical conductivity 
where C is a constant depending upon geometry. The relationship between 
attenuation factor and the shield time constant is simple for a shield 
whose thickness is less than the magnetic skin depth for the frequency of 
the asynchronous magnetic fields. If the shield is thicker than a skin 
depth, the expression is less direct. 
For (wT, )? > 1, o. = angular frequency of asynchronous flux, 


relative to the rotor, 





(II-1) 


2 \% 
6 «= Os skin depth at frequency ws 


A = shield thickness 


R,R, are the inner and outer shield radii. 


The above expressions do not correspond for 6 = A because higher- 
order correction terms have been omitted from the second expression for 


6 <A. However, it should be clear that, as the shield thickness A or 
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the shield conductivity O increases, TS, increases, and the attenuation 
factor becomes a smaller number. 

Since protection of the field winding against conditions tending to 
drive it normal is related to the attenuation factor, the first impression 
is that Katt should be made as small as possible. This implies a shield 
with a very long A However, the shield must also serve as a damper 
winding. A shield with a long time constant has a low resistance, and is 
not very effective in damping machine swings. Reference [5] has shown 
that the optimum value of shield time constant for damping is given by 


x +x 


w T = 


a's," x"¥x, a 


where Ww. is the angular frequency of machine swings, and Xx, is the exter- 
nal reactance through which the machine is connected to an infinite 
electrical bus. Reference [5] shows that this optimum shield time cons- 
tant is nearly independent of Xx, because of the dependence of W, upon x,- 
Therefore, given the synchronous and subtransient reactances, the shield 
time constant for optimum damping is fairly well determined. As eet is 
increased or decreased from this optimum value, the damping coefficient 
is decreased, as predicted by reference [5]. 

Shield time constant and attenuation factor are related, as was 
indicated previously. Therefore the shield time constant cannot be set 
arbitrarily at the optimum damping value without consideration of the 
requirements for attenuation. These will be set by the induced field 
current and thermal dissipation in the cryogenic region during faults 


which can drive the superconductor normal. To provide a basis for the 


tradeoff consideration between damping and attenuation, induced field 
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current and heating are covered in the following section. 

Induced Current in the Field Winding 

Type II superconductors carrying direct currents with no loss will 
exhibit a transition to normal conductivity if sufficient alternating cur- 
rent is superimposed on the direct current. Losses are known to be 
associated with the alternating currents, [11,12], but evidence does exist 
that the transition is not independent of the de level; that is, the 
transition appears to occur if the peak current (ac and dc) attains a 
critical valve [10]. 

Following any type of short circuit on a synchronous generator, 
there is a rise in the direct current level in the field winding, and a 
superimposed alternating current. These changes in field current are due 
to the demagnetizing effect of the armature fault currents and the mag- 
netic flux trapped by the shorted armature phases. 

A conducting shield surrounding the field winding will attenuate the 
level of alternating current induced, but can only limit the rate of rise 
of the direct current level. In section IV-4, expressions for the induced 
field currents, ac and dc, for three-phase and line-to-line short circuits 
on a generator are presented. Following a three-phase short circuit on 
a generator occurring at t = 0, the field current will have the following 
form, assuming that the field time constant is much longer than the shield 


time constant: 
t 


x- x' ~ eu 
i iit +e (fl -e d )j +i 


£ f£ 


ac 


The field current is expressed in per-unit, where the base is the 


field current necessary to produce rated open-circuit voltage. The 
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value i, is the per-unit field current before the fault; x,x' are the 


f 


d is the short circuit 


subtransient time constant which is related to To by Eq. (A-15); i is 


fac 


the peak per-unit value of induced alternating current in the field wind- 


synchronous and transient reactances in per unit; T 


ing after the fault. 

The field current before the short circuit Leo will increase by an 
amount Ai, = (x-x" ALS due to the demagnetizing effect of the armature 
currents. (For the three=phase fault Ai; = V,,/x', but other cases can 
be evaluated using Table IV-1.) This rise occurs exponentially with the 
subtransient time constant (which is related to the shell time constant). 
The maximum level reached by the field current will depend upon the time 
necessary for the circuit breaker to interrupt the fault, and upon the 
value of the shield time constant. If we assume that the circuit breaker 
can operate in 5 cycles, the maximum per-unit increase in field current 
can be calculated as a function of the shield time cmstant. For a sub- 
transient time constant short compared to 5 cycles, the field current rise 
will be about 0.6 per unit for a three-phase short circuit upon the 1000 
MVA machine of section II-4. Other cases can be slightly worse. For 
example, operation at zero power factor overexcited with a line-to-line 
short circuit produces a field current of .84 per unit. Other cases can 


be estimated, using Table IV-1. Unsymmetrical short circuits from load 


depend on initial conditions and system reactance. 
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The point of the last paragraph is that, for a representative fault 
condition, such as a three-phase short circuit, it is possible to relate 
field current rise to shield time constant. The longer the shield time 
constant, the less the rise. But in making the shield time constant 
longer, damping ability of the shield is degraded. From reference [5] it 
is possible to calculate the change in damping coefficient for an increase 
in shield time constant above the optimum value. In addition, the level 
of induced alternating field current can be calculated, using Eqs. (III-8), 
Table IV-2, and Eqs. (A-53-54). 

The results are plotted in Fig. II-l. The ratio Kp/Kp ant relates 
damping coefficient to its optimum value and the ratio Tyl ten ope relate 
shield time constant to the value of To producing optimum damping. The 


plot of Kp/& ane versus T/T indicates the loss of damping versus 


S1 opt 
increase (or decrease) of shield time constant from its value for opti- 
mum damping. Also plotted are the per-unit values of field current rise 


and peak alternating induced field current versus T_ /T for a three- 


S S81 opt 


phase short circuit from load with breaker opening after 5 cycles. For 


example, if T. is increased to twice its optimum value, damping will be 


1 
80% of optimum, field current rise will be limited to 39 per unit, and 
peak induced alternating field current will be 0.25 per unit. 

Figure II-l is intended to give an indication of the loss of damp- 
ing associated with limiting the field current rise during a fault. The 
requirements of machine damping will depend on the installation. The 
field current rise will determine how the operating point of the super- 


conductor is set so that superconductivity will not be lost before the 


breaker opens. It is obvious from Fig. II-l that, even for the best 
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FIGURE II-1 


Damping and Field Current Surge vs. Shield Time Constant 
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damping shield, the induced ac is small and the dissipation will probably 


have a smaller effect in causing a transition than the current level rise. 


Heating in the 4.2° K Region During Faults 


During the period immediately following a fault, considerable heat- 
ing occurs within the rotor, due to trapped armature flux and imbalance 
of the shorted phases. The heating which is of most immediate interest, 
as far as maintenance of superconductivity is concerned, is the heating 
which occurs within the 4.2° K region where the field winding is located. 
There is a steady heat leak into this region, due to conduction through 
the mechanical supports from the room-temperature region radiation from 
the thermal shield, and dissipation from asynchronous flux which pene- 
trates the electrical shield. 

The cooling system must be sized to handle this steady-state energy 
input rate. During the fault period, an increase in this 4.2° K region 
heat load is experienced, the magnitude of which is dependent upon the 
electrical shield's attenuation factor. The sources of this increase are 
the larger asynchronous flux during the fault, and the increased radiation 
from the shield due to its temperature rise. 

Two cases are of interest concerning heating in the 4° K region. 
First is the short-term fault which is cleared in 5 to 15 cycles. The 
total energy input to the 4° region must not raise the fleld winding 
above its critical temperature, so that superconductivity is not lost. 
Second is the case of a sustained fault which cannot be iaterrupted by 
the circuit breaker. In this case, the primary concern is not necessarily 


maintenance of superconductivity, but prevention of permanent damage to 
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the machine. If the field winding is driven normal, a large dump of magnetic 
energy will occur within the cryogenic region. This must be anticipated 
in order that the rotor be designed to withstand such an occurance. 

The rate of energy dssipation in the containment vessel depends on 
the attenuation factor of the electrical shield, and on the thickness of 
the stainless steel vessel. For the 1000 MVA machine of section II-4a, 
stainless steel thickness of .29 inches gives maximum dissipation. A 
thinner shell is resistance-limited; a thicker one is reactance-limited 
until the shell is thicker than a skin depth, which is about 1.8 inches. 
The short term (5 cycles) temperature rise in the 4° region is shown in 
Fig. II-2 as a function of the electrical shield attenuation factor, and 
with containment vessel thickness as a parameter. This temperature rise 
depends upon the thermal capacity of the rotor. The liquid helium has 
by far the largest specific heat of any material in the 4.2° K region, 
and contributes almost all the total heat capacity. It has been assumed 
here that helium occupies 10% of the field winding volume. A larger per- 
centage would give more heat capacity; however, if the field winding goes 
normal, all the liquid will be vaporized and must be vented. Therefore, 
for safety reasons, it is probably advisable to keep the percentage of 
helium as low as possible. 

The results shown in Fig. II-2 indicate that maintenance of supercon- 
ductivity for a short circuit cleared in five cycles will require that 
the attenuation factor for frequency WO, be kept below .01 in order to 
keep the temperature rise less than one degree. Optimum damping can still 
be realized for an attenuation factor less than .01 if the shield is 


thick compared to a skin depth, i.e., when the second part of Eq. (II-1) 
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FIGURE II-2 


Temperature Rise in 4 Region vs. Shield Attenuation Factor 
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applies. In section II-4 is presented an example in which optimum damp- 
ing is obtained with an attenuation factor of .0016. 

Another interpretation of Fig. II-2 is to say that the containment 
vessel should be made of a nonconducting material, or that some provision 
be made for breaking up the long eddy current paths in the containment 
vessel. However, stainless steel appears to be the best structural mater- 
ial for the containment vessel, and design of a laminated vessel would be 
difficult. It is probably more practical to assume that the proposed con- 
figuration of a stainless steel vessel will be used. In such case, an 
electrical shield with adequate attenuation will be required. Additional 
benefits result from an attenuation factor of .01 or less, as will be 
shown in the following section concerning structural integrity of the 
field winding. This attenuation can be realized consistent with reason- 
able damping. Therefore, it remains to be demonstrated that a shield 
can be designed to produce this attenuation and damping, and have suffi- 


cient structural and thermal capacity to survive fault condition itself. 


7 


— 
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TI-2 Maintaining Structural Integrity of the Rotor During Faults 


During the period following an armature short circuit, very large 
stresses are exerted on the rotor. The parts which are subjected to 
these stresses are at the same time being heated by induced currents, with 
the result that their structural properties are degraded. Maintenance 
of structural integrity involves consideration of the thermal capacity 
of the rotor parts subjected to the stresses and heating. The following 
sections discuss the magnitudes of the stresses exerted on the field wind- 
ing and the electrical shield. The thermal capacity of the rotor for 
sustained faults is then investigated. 

Rotor Fault Torques 

The torque exerted upon the rotor following a short circuit has a 
large oscillating component and a much smaller unidirectional component. 
The maximum torque occurs during the subtransient period; if the rotor 
can survive the first few cycles, it will survive the fault, provided 
the mechanical natural frequencies of the turbine generator rotor are 
not excited by the oscillating stresses. A structurally sufficient 
model must therefore be substantial enough to resist the subtransient 
stresses and have its natural vibrational frequencies above the excita- 
tion frequencies. 

Field Winding Fault Torques 

One of the principal benefits of the electrical shield is the 
attenuation of oscillating torques applied to the field winding. Table 
II-1 shows the peak per unit total torque upon the rotor and the torque 


on the field winding for several types of faults from various conditions 
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of loading. The supporting calculations are from Appendix I with the 
analytical expressions from Chapter IV. For attenuation factors less 
than .01, as were seen to be required for maintenance of superconduct- 
ivity during short term faults, the level of alternating field winding 


torque is less than 3 & of rated torque. 


- Table II-1 Fault Torques 


Peak per unit Peak Torque upon Field Wind- 
mepcs tion Torque on Rotor ing for Various Attenuation 
Factors K.., (w) 
001 01 ol 


Three-phase short circuit 
at terminals from load 


at pp = «85 6.7 0025 025 25 


Line-line short circuit 
at terminals from rated load 
at Pe = .85 10.5 0017 O17 7] 


Closure out of phase by 


The maximum torque which will be exerted upon the field winding is 
essentially the rated torque of the machine. The benefit of field wind- 
ing fault torque limitation is that the structural supports for the field 
winding can be minimized. The conduction heat leak into the 4.2° K 
region can be reduced, thus reducing the capacity of the refrigeration 


system required. 
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Fault Torques Upon the Electrical Shield 


The total reaction torque upon the rotor can be calculated indepen- 
dent of the details of damper shielding, as is done in many references 
(e.g. [14]). The total torque acting on the rotor is exerted primarily 
on the electrical shield, if the attenuation factor is low as is required 
to maintain superconductivity during the short-term fault. The peak 
alternating torques in per-unit acting on the shield are calculated from 
the expressions in section A-l-b and are listed in Table II-l for various 
fault conditions. 

If the electrical shield is to withstand the fault torques, the 
yield stress in shear of the shield material must not be exceeded at the 
maximum temperature reached by the shield under the fault condition. 
Torsional buckling of the shield is also possible, but for practical 
shield designs is a higher limit than yielding in shear. Torsional oscil- 
lation natural frequencies of practical shields are always considerably 
above the double line frequency of the fault torque for unbalanced faults. 
In section II-4 of this chapter, a shield design is proposed which is 
sufficient to resist these modes of failure. 

Since most of the fault torques are taken by the shield, its sup- 
ports must have considerably more structural strength than the field 
winding supports. It is therefore more difficult to provide thermal 
isolation of the electrical shield consistent with adequate mechanical 
support. Although the MIT 2 MVA machine has a 20° K shield rotating in 
a vacuum relative-motion gap, it appears, based on structural support 
of the shield, that a room-temperature shield would be more practical 


for power-system-size generators. This assumes that adequate steady-state 





38 
cooling of the shield under slight imbalance is possible. This will be 
considered in the section on steady state performance. However, the most 
practical configuration for the 1000 MVA machine rotor appears to be a 
heavy, well eeaponted. room-temperature electrical shield with the 4.2° K 
region containing the field winding, supported by thermal distance pieces 
which need to support only slightly above rated torque and can minimize 
the conduction heat leak. A secondary shield can be provided at approxi- 
mately 20° K to intercept thermal radiation from the room-temperature 
shield. The conduction heat leak through the field winding supports can 
also be intercepted at 20° K, making the heat load in the 4.2° K region 
quite small. 

The amount of structure necessary to support the transient torques 
upon a 1000 MVA machine is considerable. Let us assume that the field 
winding is supported by a structure adequate to withstand rated torque, 
which is about the maximum it will feel as was shown in the last section. 
If the field winding is supported by a solid shaft of stainless steel 
with yield stress of 60,000 psi, it must be seven inches in radius. To 
allow for the helium plumbing and current leads, the support must, in 
fact, be hollow. Figure II-3 shows the required field support thickness 


A, as a function of the outer field winding support radius F,. As Fy 


f 
increases, the required A, decreases, lowering the heat leak to the 4.2° F 
region. For a vaiue of E = 9 in., the rate of gain in thermal isolation 
is decreasing. If a 2-inch buildup of superconductor is assumed, the 
inner radius of the shield support is set. The righthand portion of 


Fig. Ii-3 then indicates the required thickness of the shield support 


to withstand the maximum transient torque, which is approximately 10 per 
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unit. The resulting shield support 
electrical shield 
is roughly 5.4"' in thickness. If the 


shield support 


shield were to be run at cryogenic 
— R A stainless steel 
O 
temperature, the heat leak through the w 
, superconductor 
supports would be at least 5 KW. ae ield winding support 


Coupled with the steady-state shield 
dissipation (as discussed in the next 


section), this would require nearly Rotor cross seeeen 
a megawatt of refrigerator power. 

The alternative is to run the shield at room temperature, cooling it 
in a manner similar to the rotor cooling of conventional machines. The 
heat leak through the field winding support can be limited to about 523 
watts, and no negative sequence dissipation is experienced under steady- 
State because of the shielding of the room-temperature shield. A 
second cryogenic shield at a temperature slightly above the field winding 
can be used to intercept the thermal radiation from the room-temperature 
shield. This second shield can also be used to provide a lower effective 
attenuation factor to asynchronous flux during faults, and thus limit the 


temperature rise during short-term faults. The over-all configuration 


will be given in section II-4. 


Normal Stresses Upon the Electrical. Shield 

During the subtransient period following an armature short circuit, 
the demagnetizing flux produced by the armature fault currents is 
excluded from the electrothermal shield. A normal stress is produced 
which is nonuniform around the shield, and which tends to flatten it into 


an elliptical shape. In addition, the asynchronous fluxes as seen on the 
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rotor produced by trapped armature fluxes and armature imbalance are 
also excluded from the shield. These fluxes also produce a traveling 
wave of deflection on the shield, so it must have sufficient rigidity to 
this type of stress to limit the maximum deflection to a tolerable level. 
If the shield is not sufficiently resistant to radial deflection itself, 
provision must be made for radial stops upon the surface of the contain- 
ment vessel inside, to limit the extent of the deflection. The ring 
natural frequencies of the shield must also be kept above the highest 
frequency of excitation, which is twice line frequency. 

The maximum radial deflection of the shield and its ring natural 
frequencies are primarily functions only of its thickness. Additional 
radial stiffness can be obtained by the addition of a stainless steel 
cylindrical shell inside the electrothermal shield. As will be seen in 
the design in section II-4, this is probably necessary to withstand the 
transient torques. In addition, radial stops located on the helium con- 
tainment vessel and supported radially through the field winding support 
can also be provided. Contact would not be made between the shield and 
the stops under normal conditions, so there would be no radial heat con- 
duction. Only during a transient deflection would the shield come into 
contact with the stops. 

Figure ILI-4 shows a plot of the maximum radial deflection under 
transient and the lowest natural ring frequency as a function of shield 
thickness. For a self-supported shield, considerable material thickness 
is required to limit resulting deflections. 


Thermal Capacity of the Electrical Shield 


In section II-l, the thermal capacity of the 4.2° K region was con- 
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sidered to determine the temperature rise for a short-term fault in order 
to maintain superconductivity in the field winding. For a sustained 
fault upon a machine which cannot be interrupted, such as a three-phase 
short circuit at the machine terminals, the primary concern is the 
ability of the electrical shield to absorb the resultant ‘heating without 
suffering permanent damage. The measure of a machine's ability to absorb 
transient heating which greatly exceeds the cooling rated is expressed 
by the quantity I *t, where I, is the per-unit negative sequence armature 
current, and t is the time for which I, is applied. For example, a 
machine with I *t = 8 seconds would withstand rT = 1 for a period of 8 
seconds without permanent damage. The industry standard for turbine gen- 
erators is I *t 2 8 seconds. The standard is presently under discussion. 
Turbine generator manufacturers are trying to have the standard lowered, 
because of the difficulty of obtaining thermal capacity in conventional 
generator designs. To satisfy this requirement, the electrical shield 
must have sufficient heat capacity to absorb energy equivalent to Tere 2 8 
seconds without damage. A value of ee ereater than 8 seconds would be 
very desirable. 

If the material and thickness of an electrical shield is selected, 
a value for it can be determined by assuming an initial operating temp- 
erature for the shield and performing a step-by-step calculation of the 
temperature as a function of time for I, = 1, taking into account the 
dependence of conductivity and specific heat upon temperature. The 
maximum temperature is determined by the loss of yield strength with 
increased temperature. At cryogenic temperature, specific heat and elec- 


trical resistivity of materials are both low; as temperature is increased, 
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both increase, but specific heat levels off near room temperature, while 
resistivity continues to increase. The result is that the Tega: associated 
with a temperature rise from 20° K to 300° K is about the same as that for 
a rise from 300° K to 650° K. 

The nature of an Bat number like 10 seconds is misleading. It 
implies that a machine can withstand a fault condition for 10 seconds, 
in which case thermal diffusion would permit heat to be distributed to 
all structure in contact with the electrical shield, and additional ther- 
mal capacity would be available from the stainless steel reinforcement 
backing the shield. In fact, what we are trying to provide is sufficient 
thermal capacity such that a terminal three-phase short circuit will not 
heat the shield to a point where structural damage can result. For a 
three-phase fault, the value of I, is of the order of (x") which is 
6.25 per unit, and the duration of the heating is approximately the arma- 
ture time constant, which is about 0.3 second. In a period of .3 seconds, 
almost no heat is conducted into the underlying stainless steel because 
of its lower thermal conductivity. The value of I *t required to absorb 


the dissipation of a three-phase terminal short circuit is 


V 
oc " 


2 
— 3 @short circuit x 


This number is at = 12.7 seconds for the 1000 MVA machine of section 
Ii-4. In order to provide this thermal capacity, only the heat capacity 
of the conducting material of the shield can be relied upon. The impli- 
cation is then that the conducting material must be made thick enough to 
provide the thermal capacity. However, from Eq. (A-56), it is seen that 
the shield time constant increases as the thickness. Calculation shows 


that a shield with sufficient thermal capacity at room temperature 
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initially, and of high-conductivity copper to reduce steady-state loss, 
would have a time constant three times the value for optimum damping, 

thus reducing the damping to 60% of optimum. For a copper shield at 20° K 
if sufficient thickness for heat capacity is provided, the electrothermal 
shield time constant will be nearly 150 times the optimum value, and pro- 
duce almost no damping. 

In order to overcome this problem, there are two possible alterna- 
tives. We can use a material of lower conductivity, permitting a thicker 
shield without an excessively long time constant. In order that the shield 
attenuate the asynchronous fields sufficiently to limit dissipation in the 
stainless steel support structure below it, the shield material should be 
thicker than the magnetic skin depth at the asynchronous frequency. The 
steady-state power dissipation is proportional to the square root of 
the shield resistivity. Hence, the higher-resistance material will re- 
sult in larger steady-state shield dissipation. 

A second alternative for providing thermal capacity of the shield, 
and limiting the shield time constant, is the use of a circumferentially 
slotted shield. This technique permits the desipner to use a high- 
conductivity material for the shield, like electrolytic copper, and reduce 
the steady-state losses. The disadvantage is that the copper has low 
yield strength and must be held down upon its supporting tube with steel 
banding. The technique of slotting does permit independent selection 
of thermal capacity and time constant with the use of high-conductivity, 
low-loss material. The slotting technique is explained in the following 
paragraphs. 


If narrow circumferential slots are machined into the shield, the 
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time constant of the shield is reduced because the resistance of the 
shield to axial current flow has been increased. In chapter V are pre- 
sented expressions for the resistance increase associated with such an 
operation, along with experimental verification. The result is that, by 
proper selection of depth and spacing of slots, it is possible to adjust 
the time constant (and hence the attenuation factor) of a conducting 
shield without seriously affecting its thermal capacity. 

Figure II-5 shows the cross- 
section of a shield taken by a 
plane containing the cylindrical 


axis. The shield's resistance - 





to current flow and hence its 


time constant, are determined by 


Slotted shield cross section 


the quantities S, t, and h, 
while the thermal capacity is Ebeuneglle? 
determined primarily by the thickness, w. 
For example, a room-temperature copper shield 1.5 inches thick in the 
1000 MVA machine would have an a: = 20 seconds, but a time constant of 


T -384 seconds. Slots .1 inch wide and one inch deep, spaced one inch 


$1 
apart, will reduce the time constant to TS, = .189 seconds, but maintain 
at = 18.8 seconds. The thermal diffusion time associated with con- 
duction of heat through the one-inch-thick "fins" is .52 seconds, indi- 
cating that, for a .327 sec. armature time constant, approximately full 
advantage of the copper's thermal capacity will be realized. However, 


the thermal diffusion time for one inch of stainless steel is 14.5 sec., 


indicating that it will provide very little thermal capacity for transients. 





&7 

Determination of the thermal capacity, or It, of a shield, depends 
on the maximum temperature to which the shield will be permitted to rise. 
For the example quoted in the last paragraph, the copper shield was per- 
mitted to rise to 650° K or 377° C. This is consistent with the tempera- 
ture rise estimated for the amortissem of conventional generators for an 
I *t of 8 seconds [20]. The rotor surface temperature for conventional 
generators for It = 8 seconds is estimated at 460° C and for r*t = 12 
seconds, at 580° C from reference [20]. 

A temperature of 377° C is in the low annealing range for copper; 
hence, the maximum yield strength which can be counted upon is about 2000 
psi for fully annealed copper. As shown in section II-4, this is insuf- 
ficient to withstand the centrifugal stresses; sufficient strength to 
withstand the fault stresses can be provided by stainless steel backing. 


The big problem is in holding the copper on with steel banding. 


Use of the slotted electrical shield proposed above can permit the 
machine designer to determine the thermal capacity, and hence the I *t, 
of the generator consistent with damping and attenuation requirements. 

It is possible to provide larger thermal capacity with the thick copper 
shield than with the fron rotor of a conventional machine, because of 
the larger thermal conductivity of copper. The tron rotor of the con- 
ventional machine is heated to extremely high temperatures before much 
heat has been conducted into the bulk of the rotor. 

The first alternative, of a thick shield of higher resistance mater- 
ial, has been chosen for the design presented in section II-4. The 


material selected for the shield (Phosphor Bronze) has considerably 


higher yield strength than electrolytic copper; it can therefore withstand 
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the centrifugal stresses without external support rings, and can carry a 
significant portion of the shear stresses resulting from fault torques. 
The steady state dissipation for armature imbalance is higher than it 
would be for a slotted high-conductivity shield, but it is still small 
compared to the rotor heat load of a conventional machine of comparable 
size. 
II-3 Steady-State Performance of the Electrical Shield for Unbalanced 
Load 

The primary design consideration concerning steady-state operation 
of the electrical shield is the power dissipation due to armature current 
imbalance. To keep this dissipation as low as possible for a given mater- 
ial, the shield must be thicker than the skin depth of the induced cur- 
rents due to imbalance, which is twice line frequency. The dissipation 
is then proportional to the square root of shield resistivity. This 
implies that advantage is gained from a cryogenic shield; however, the 
important factor is the amount of refrigerator power required to remove 
the heat dissipated in the shield. Added to this negative-sequence heat- 
ing is the conduction heat leak for a refrigerated shield. Because of 
the heavy supports necessary to support transient torques, this heat 
leak is quite large. 

Table II-2 gives a comparison between the heat loads due to negative 
sequence and conduction for shields at room temperature, 77° K, and 20° K. 
Also computed are the refrigerator power for the cryogenic shields and 
the cooling power for the room-temperature shield. The numbers concern- 
ing refrigerator power should be taken only as qualitative, but do indi- 


cate the very definite trend that higher shield temperature can be cooled 
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more economically. 


Table I[I-2 


Cooling Power Required for Shield as a Function of 
Shield Temperature 


Shield Temp- Conduction Dissipation Total : Total blower 
erature Heat Load for I = .1 Heat Load refrig. or compressor 
°K (KW) (KW) (KW) Power (KW) 
20 3 10 ‘5 50 750 
77 4 35.4 39.4 9.4 370 
300 0 150 150 (add 100KW 250 


for fan power) 


The other steady-state considerations of the shield's performance 
are its mechanical natural frequencies, and the centrifugal stress level 
in the conducting layer. The critical natural frequencies include the 
torsional frequency of the shield, the ring natural frequency for radial 
vibration, and the lateral vibration frequency of the shield between 
its end supports. The formulae for these calculations are given in 
chapter III, and results for a 1000 MVA machines are in section II-4. 

None of these natural frequencies appears to be a problem because of 
the thick shield support required by transient torques. 

The stress in an outer, high-conductivity copper layer due to centri- 
fugal stress is a problem, because of the rotor diameter. The yield 
stress of the copper is low because of its purity and annealed condition. 
The critical radius for centrifugal stress at a yield stress of 10,000 psi 
is 13 inches. Since the rotor must be larger than this, it would be nec- 


essary to provide retaining rings of stainless steel upon the copper 
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shield. The dissipation in these rings is small because of their 
limited axial length. Another possibility is to use fiberglass or some 
other nonconductor to tie down the electrothermal shield copper layer. 
Use of a lower-conductivity, higher-strength material can eliminate this 


problem, as has been done in the design presented in section II-4. 


II-4 1000 MVA Superconducting Generator; an Example 
Designs for 1000 MVA superconducting generators by Thullen [17] and 


Kirtley {21] have not included consideration of the electromechanical 
stresses imposed on the rotor during faults. Least-weight designs have 
been developed, but calculations show that they are not structurally suf- 
ficient in the rotor to withstand the transient torques and forces. The 
design presented in this section has started with the least-weight con- 
figuration developed by Kirtley [21] and altered it in such a manner as 
was necessary to produce a machine which can survive fault conditions 
consistent with the results of this study. 

Considerable weight has been added to the rotor; however, this is 
necessary to provide structural integrity. The least-weight machine 
on the basis of Kirtley's work, had been shown to be long and small in 
diameter, as opposed to the previous design developed by Thullen [17]. 
It has been necessary in the present design to increase rotor diameter 
somewhat; however, an effort was made to keep the rotor diameter as 
small as possible, in the belief that this is probably nearer the least- 
weight configuration that could be developed consistent with fault- 
survival ability. No weight optimization has been done for the present 
design, but it can be done using Kirtley's techniques with the requirement 


for fault survival included. The design used by Kirtley did not have an 
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iron shield over the end turns. The end turns are then assumed to con- 
tribute inductance, but generate no voltage. More recent designs have 
included iron shielded end turns which contribute more to effective 
length. Therefore, optimum designs may in fact be shorter and larger 

in diameter. Such a machine makes satisfaction of fault-survival criteria 
even easier, because of the increased effectiveness of rotor structure 

at a larger radius. 

It was decided to develop a design which maintained a subtransient 
reactance of .15 per unit. Decreasing the subtransient reactance to a 
lower value results in a rapid increase in the maximum per-unit transient 
torque which must be tolerated. A 152% subtransient reactance results in 
a maximum subtransient torque of about 10 per unit; this is for a line- 
to-line short circuit at the terminals while operating at rated condition. 
The three-phase short circuit at the terminals results in a maximum 
torque of 6.7 per unit. The case of 10 per unit torque may seem rather 
drastic, but the object is to show that it is possible to develop a design 
consistent with such a requirement; the occurrence of such a fault is 
entirely conceivable. 

The thermal capacity of the shield for this design will have a 
thermal capacity indicated by an It of 20 seconds for a 361° C temper- 
ature rise of the shield. This is considerably more than the available 
oe for conventional gnerators of similar rating. 

Maintenance of superconductivity through the short-term fault de- 
pends upon the shield attenuation factor and the thermal capacity of the 
liquid helium region. The attenuation resulting from the electrical 


shield for this design is .0016; the temperature rise in the helium region 
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for a short circuit cleared in five cycles less than .1° K, and for one 
cleared in 15 cycles is less than .3° K. The induced alternating field 
current is .0048 per unit, and the dc rise is .24 per unit. Proper place- 
ment of the superconducting winding operating point relative to short 
sample characteristics should guarantee maintenance of superconductivity 
through these conditions. The secondary, or thermal, shield can be de- 
signed to decrease the problem of field current surge. Since this shield 
is subjected to almost no electromechanical stresses, it can be designed 
with the proper time constant necessary to limit the field current surge 
consistent with the ability of the exciter to control the field arrent 
change. 

The damping coefficient of the present design has been optimized for 
the estimated swing frequency of the turbine-generator system. The amount 
of damping will depend on the shield temperature during the period of 
machine swinging. However, the power dissipated in swinging will be dis- 
tributed throughout both the shield and its stainless steel support tube. 
Considerable thermal capacity is therefore available, and the temperature 
rise should be rather small. Therefore the effective time constant and 
damping coefficient should not change radically throughout the period 
of machine oscillations. 

Mechanical natural frequencies of the rotor relative to lateral 
vibrations, torsional oscillation, and ring vibrations are all high 
enough that they should cause no problems. The field winding supports 
have been minimized, making the lateral natural frequency of the field 
winding the lowest expected natural frequency. This calculation has 


neglected any stiffening due to the winding itself; the true lateral 
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natural frequency will therefore be somewhat higher. Steady-state 
dissipation due to 10% negative sequence is 141 KW. The windage will be 
comparable to a conventional machine. However, the shield dissipation is 
less than 5.2% of the field winding iR losses for a conventional generator 
of this size. Therefore, cooling of the shield by conventional techniques 
should not be difficult; the rotor could operate in a partial vacuum to 
eliminate windage, but shaft seals must then be provided. 

The heat load in the 4.2° K region under stady-state conditions is 
primarily conduction through the supports and radiation from the secondary 
shield. These losses can be intercepted at a slightly higher temperature, 
such as 20° K, reducing the 4.2° K heat load to a few watts. The second- 
ary shield will act primarily as a radiation shield, and will not be 
subjected to the subtransient mechanical stresses and heating during 
faults whicn is taken by the primary, room-temperature shield. The cryo- 
genic secondary shield can be made quite thin, of high-conductivity mater- 
ial, and cooled by conduction from the ends. 

Figure II-6 snows a drawing of the proposed design, and Fig. II-7 
shows more detail of the rotor. Table II-3 lists machine parameters and 


characteristics. 
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FIGURE II-~7 


1000 MVA Design Rotor Cross Section 
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Table II-3 


1000 MVA Superconducting Generator 


Stator Dimensions: 


Inner Radius of Conductors 

Outer Radius of Conductors 

Inner Radius of Iron Magnetic Shield 
Outer Radius of Iron Magnetic Shield 
Length Straight Section 

Length Each End-turn Region 


Rotor Dimensions: 
Inner Field Support Radius 
Outer Field Support Radius 
Inner Superconductor Radius 
Outer Superconductor Radius 
Helium Containment Vessel Radius 
Secondary Shield Radius 
Inner Primary Shield Support 
Outer Primary Shield Support 
Inner Primary Shield Radius 
Outer Primary Shield Radius 


Reactances: 
Synchronous Reactance x 
Transient Reactance x' 


Subtransient Reactance x" 


1359 
25.0 
29°0 
41.7 
121.3 
22.5 


8.95 
10.0 
10.0 
12.0 
12.0 
D2 
12.2 
16.0 
16.0 
17.6 


237 
215 


in. 
oleic 
in. 
in. 
in. 


in. 


in. 
in. 
in. 
in. 
in. 
in. 
in. 
in. 
in. 


in. 
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Primary Shield and Support 
Shield Material: 
Phosphor Bronze A 
954 Cu, 54 Sn 
Hardness: Rockwell B-91 
oy % 80,000 psi ; Fatigue Strength for 10° cycles = 33,000 psi 


90 (300° K) = .107 ; 6 (w,) = .835", 5 (2) = ,59" 


Shield Support Material: 
Stainless Steel ie ~ 60,000 psi 


p(300° K) = .725 ; 6() e217" , “ecco, #*1.54" 


Thickness Time Constant (ms) 
Shield Primary ik6"" 159 
Support Material 5.5 55.5 


214.5 ms. 


Ds optimum damping = 212 ms for HW. = Oe: le, f. = 1.6 hz. 


Power Dissipation in Phosphor Bronze: 


rT = .1 141 KW 


K (w_) af 
att Oe .0207 Kao 65 .076 
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Stresses in Shield 


Syield stainless steel = 60,000 psi 


Syield phosphor bronze = 80,000 psi 
Max. Transient Torque = 10 = 240 x 10° in. bf. 
Oo Zeal R 
shear due to transient torque ; oO = 3eNMeuo psi 
mR a R;) 


o, (centrifugal hoop stress) = 18,340 psi 


Max. shear from Mohr's circle 


Se max = 37,500 psi = .47 


tensile yield 
Overspeed centrifugal stress (5,400 rpm) 
o (5400 rpm) = 41,100 psi 
t 
Cig 2a 20,600 psi 


Maximum Radial Deflection for 3¢ Short Circuit = .037 in. 
Mechanical Natural Frequencies 
Lateral Vibration of Shield 1 ls 176 hz. 
Lateral vibration of Field Winding £ = 50 hz. 
Torsional Natural Frequency Shield ES oe 442 hz. 


Torsional Natural Frequency Field 
Winding f, = 256 hz. 


Ring Frequency at Standstill Shield = 655 hz. 


Ring Frequency at Standstill Field 
Winding eS 








a 


Thermal Capacity 
Shield Initial Temperature = 27° C 


. a Phosphor Bronze 
I, t max 
10 sec. 193° C 
20 sec. 361° C 
25 sec. 446° C 


Max Shield Temp. for three-phase short circuit 
at terminals = 251° C 
Max Shield Temp. for line-to-line short circuit at terminals 


assuming field winding discharged in 5 sec. = 273° C 


Heat Loads in Cryogenic Region 


Secondary Shield at 20° K 


Thermal distance piece length = 16 in. 
thickness = 1.05 in. 


Conduction 523 watts 
Radiation 103 watts 
Negative sequence ( < 1 watt) 0 


626 watts @ 20° K 
Field Winding Region 4.2° K 


Thermal distance piece lenpth = 8 in. 


thickness = 1.05 in. 


Conduction 2 watts 
Radiation neg. 
Refrigerator Compressor Power 33.3 KW ( < 10° of machine 


rating) 





60 


Estimated Weights 


Rotor: 

Field winding support 320 

Field winding 5070 

Shield Support 14680 

Shield 9190 

Shaft (20 inches each end) 8060 

37320 lbs. rotor 

Stator: 

Armature 24970 

Iron Shield 136700 


161670 lbs. stator 


198990 lbs. total 


Rating 985 MVA -202 lbs/KVA 
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The proposed design has a synchronous reactance of .59 com- 
pared to nearly 2.0 for a conventional machine of similar size. The .202 
lbs/KVA weight is nearly an order of magnitude less than that for a con- 
ventional machine. The thermal capacity, as indicated by ee is at 
least 20 seconds compared to 8 to 10 seconds of present iron rotor design. 
The efficiency is about 99.8% compared to 98.7%, and the rotor natural 
frequencies are all high compared to conventional machines. Thus, the 
advantages associated with superconducting generators have been main- 


tained while satisfying the criteria for survival of fault conditions. 
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CHAPTER III 


Steady-State Effects Upon the Electrical and Thermal Shields 


1) Heating of the Rotor 


A heat load is imposed upon the rotor of the generator by eddy currents, 
conduction through structure, and radiation. The heat load will be distri- 
buted among the electrical shield, the cryopenic thermal shield, the stain- 
less steel can and supports, and the superconducting winding, depending on 
the design of the electrical and thermal shields. In chapters III and IV, 
the word "electrothermal shield'' is sometimes used. This refers, in gen- 
eral, to the electrical and thermal shield combination. When considering 
mechanical stresses, attenuation, and eddy current heating, it refers to 
the electrical shield; for thermal radiation, it refers to the thermal 
(cryogenic) shield. 

a) Sources of Heating 
1. Radiation 
The heat flux due to radiation between two coaxial cylinders 


assuming diffuse reflection, is 


P| 
i 
Q 


(T.* = i ) (III-1) 


ie 


- 1) 


m 
+ 
Eo 
on 
ie 


> 


where ae are the emissivities of surfaces having areas A, a A, 


and temperatures TY + T,- The areas A, + A, are nearly equal, and 


Te ss ee 


2 
c. re ~E CE = area of thermal shield 


1 2 


T, = stator temperature ; oO = 5.33 x 10 ° watts/ft? °K 
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(2) Stator Phase Imbalance 

The load supplied by most power system generators is not perfectly 
symmetric with respect to the three phases resulting in unbalanced cur- 
rent between the three phases. The imbalance of currents results in 
a negative sequence component of armature current, that is, one producing 
a magnetic field which rotates in the reverse direction to the rotor 
rotation. Formulae for determining the magnitude of the negative sequence 
current component are well known, and given in any power systems text- 
book. The newest standards for power system generators require that 
they be able to withstand 10% negative sequence, steady state. 

The power dissipation in the shield due to a per-unit negative 
sequence current I, is obtained from Eq. (A-50) and Eq. (A-52), where 


W = 20 andn#l. 
oO 


Negative sequence Power Dissipation: 


R2 (2 wT) 
a AT BTS (See fe 8s! 
fi @*) (=) "| Oo Si 


a (III-2) 
chi rey Br Kare 
(obey 
Oo oO 
BL, * “o°e {l-x + i a-x)(52) 3 ee 


J, = rated armature current density amp/m* rms 
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(3) Stator Harmonics 
If all three phases of the peaoaator are perfectly sinusoidally 
wound, a uniform magnetic field will be produced rotating at synchronous 
angular frequency by balanced phase currents. Spatial harmonics of 
order higher than the fundamental produce magnetic fields which rotate 
at different angular frequencies, some of which will induce currents in 
the shield and produce losses. 

Harmonics which are multiples of three produce no rotating component; 
harmonics 5, ll, 1/7, .... 6m+} produce negative sequence fields rotating 
at angular frequency w/n, or at w +=) relative to the rotor; har- 
monics 7, 13, 19, ... 6m + 1 produce positive sequence fields at angular 
frequency w @ - +) relative to the rotor; since the windings are sym- 
metric odd distributions about their axes, there are no even harmonics. 

The harmonics for a distributed winding between radii Sy and 5, 


with 60° phase winding angle are 


J = = ee S 2nt 
n nT co 3 
Gy Tt )* 
ae = of 2B nace (ILI-~3) 
+ Gy "y ne ci + (ote n 
R 2 
“it t & ay “on oe 


, . 2 4 
mo” 


The relative angular velocity of the nth harmonic Won is given by 
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1 
W = w t =) 


ER oe ee 7 om 1 


2 
‘ 


jee ees OT L 


1 
oC _ n 


n+2 


Ho Roe 2—n 2—n Ce a n+2 
=|— = + _ ¥ 
De 2°| 3 =n 56 (1 x ) = 728 ne (1 x ) J (IIE 4) 











The armature is usually made up of discrete bars which yield large 
amplitude spatial harmonics for n near the number of bars per row. 
Some concern that losses due to these bar harmonics might be considerable 


led to the following analysis. 








Each row of armature bars can be modeled 
as impulses of current at a constant radius. 
The magnitude of the nth harmonic current 
sheet representing the row of bars shown 


is calculated from the formula 


g 

: 6 ¥21 

——— 2 sinn@,;  (III-5) 
a | 


I = current in each bar rms 
2 = one half the number of bars on one 
side of the machine in the row at r= S. 


Uo Ram! Sue... 
Be x ©} a +(F) 1K, (III-6) 


The power dissipation due to each harmonic for each row of bars can be 


evaluated by Eq.s. (III-6) and (III-3) using the appropriate value of n 
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and S for each harmonic and row. The total time average power is the 
sum of the contributions of all harmonics for all rows, because the har- 
monics are orthogonal and the same harmonic for different rows are in 
phase. 
b) Distribution of the Rotor Heat Load 

The rotor of a superconducting generator will be arranged with the 
Superconductor mounted upon a center spool, probably of stainless steel, 
a stainless steel containment vessel around the field winding to contain 
the liquid helium, with a thermal shield at approximately 20° K and an 
electrical shield at room temperature. If no electrical or thermal shield . 
were included, the steady-state negative sequence component would induce 
currents in the stainless steel containment vessel and impose a heat load 
in the 4° K region. Since the stainless steel shell is a poor shield, 
due to its poor electrical conductivity, alternating currents would also 
be induced in the superconducting winding, producing loss there. The 
amount of dissipation within the superconducting wire and its consequences 
can be determined only with some operating experience with the experimen- 
tal machine being built. However, the magnitude of the losses in the 
containment shell can be predicted, based on the electrical shield para- 
meters, and the magnitude of the induced field current can be predicted. 

1) Electrical Shield Heat Load 
The total electrical shield heat load will be composed of 

the sum of negative sequence heating and stator harmonic heating. The 
calculation of these components is made using Eq. (III-2) for negative 
sequence and (III-3) for armature harmonics. The calculation of the 


resultant steady-state shield temperature must involve the cooling 
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mechanism for the shield. For example, the M.I.T. 2 MVA machine has a 
cryogenic shield cooled by conduction at the ends. Some heat will be 
lost by radiation to the colder region inside. This component of cooling 
is really more important in predicting the heat load on the radiation 
shield inside; it can be calculated from Eq. (III-l) after the approxi- 


mate operating temperature of the electrical shield has been determined. 


2) Stainless Steel Containment Vessel Heat Load 
Radiant energy from the thermal shield and asynchronous mag- 
netic fields which penetrate the electrical shield impose the heat load 
on the stainless steel vessel. If the shield is effective, none of these 
loads should be significant. The negative sequence heating will be the 
largest. The negative sequence magnetic field seen at the stainless 
steel vessel will be the negative sequence field seen by the electrical 


shield reduced by the attenuation factors presented in Appendix I. 


"12s Wer “att 

Bios = flux density at stainless steel shell (assuming 
the shell is resistance-limited) 

Bio = from Eq. (III-2) 

a = from Eq. (A-53) or Eq. (A-54) 


The time-average power dissipated in the stainless steel shell is 


Summarized as 
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P = 47x? w *ho Q Be 
Sewsinless steel - e° 12s 
To 1 as ee 2) 
Bios ene; (1-x) GP een are 
(III-7) 
= 2 
2 A <6 
Key = fugl + (2 T, ) 
2v72 2 ee A/6 A>6é 
x = radius of stainless steel shell 
h = thickness of ut " 
Ohi. = electrical conductivity of stainless steel shell 
i = time constant of electrothermal shield 
A = thickness of electrothermal shield 
6 = skin depth in electrothermal shield at frequency Zu). 
l 
- (54 \* 
0 Ho 
o = conductivity of electrothermal shield 


Equation (III-7) appears to neglect the fact that the radial field at 
the stainless steel shell is reduced by the currents induced in the 
field winding. The flux produced by the ac field currents is oscillat- 
ing, and can be resolved into two counter-rotating components relative 
to the stainless steel containment vessel. One of these components 
reduces the imposed negative sequence field, and hence reduces the 
losses; the other, of equal magnitude, rotates in the opposite direction 
relative to the shell, producing additional losses which are equal to 


the amount by which the losses due to the imposed field were reduced. 
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Hence, the total stainless steel vessel losses are as given by Eq. (III-7) 


even when the effect of the field winding is included. 


III-(2) Alternating Currents in the Field Winding Due to 


Negative sequence 


The magnitude of the ac induced-field current depends on the ac 


direct-axis armature current, according to Eq. (A-18) 


M 
df 1 
pl pe ean (I. ) 
tt 
"io FEF a (Tox? ah 








ac 


where Ww is the frequency of the direct-axis ac current. For a steady- 
state negative sequence current, the direct-axis ac current equals 

the magnitude of the negative sequence current. This equation was 
derived for a single-time constant damper winding. The field current 
expression can be written in terms of the shell attenuation in order 

to generalize for thick as well as thin shells. The attenuation factor 
is evaluated at frequency 2° 


For resistance-limited shells, (2 T55%< 1), 


= (x, - x') K I 


|| d’? “att ~2 
ac 
For reactance-limited shells (2u Lio >> 1) (III-8) 
eee = ae 
edd 
| |. AEE? 2 


(xi 7 x4) 


The expressions for Katt are given in Eq. (III-7); all quantities in 


t 
Eq. (III-8) are in per unit. 
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III-(3) Steady-State Mechanical Excitation of the Electrical Shield 


The magnetic field produced by steady-state negative sequence cur- 
rents will be excluded from the region inside the E.T. shield. The 
interaction of the imposed field and the induced shield currents will pro- 
duce a traveling wave of normal force moving relative to the shield. 

The shield should be designed so that its ring frequencies as derived 
in Sec. A-3 are not near the driving frequency of these normal forces. 
The normal electrical force per unit area is given in section A-2: 


eoecesndee He inside] 


o = i ulH -H . 
- = Neato ~~ a (ra) 


“A Uo 

17 (Me ~ ide 

r outside inside 

If the shield is effective in intercepting the negative sequence magnetic 


maeld, H will be much less than H . This is equivalent to 
6 
inside outside 


requiring that OT >> 1 where Ww. is the angular frequency of the travel- 


ing field with respect to the rotor. For negative sequence fields, 


£f He m is neglected, the force can be written 
side 


Oo i. 9 U1 He| cos (wt - 6) 
(ILI-10) 


S Sa 
oO 1 3 oO 3 Te 


At a fixed point on the shield in the rotor frame, the alternating com- 


ponent of force will vary with angular frequency 2 Ws For negative 
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sequence W_ = ZW, and the force will vary at frequency 4 Wo Therefore, 
the lowest ring natural frequency should be above 4 We 

The maximum deflection of the shell for a constant load of the form 
F cos? §@ is calculated in Eq. (A-66). The steady-state deflection for a 
mechanical system with natural frequency W and no damping, driven at 


frequency w, by a force P, sin w,t will be 


d d 


1 @ 
x = xX cr Sin Wt 
Sine Wa 
ay 


where Xo is the deflection for a constant force Pye An upper bound on 
the steady-state deflection of the shell, due to imbalance, is obtained 
by assuming no damping of the ring vibrational motion. The radial deflec- 


tion of thé shell for negative sequence current I, is given by Eqs. 








Mii-1 1): 
1 2 
Oe tse 
372 3 S S 
= eee OMe. ty t 1-3) {22 
to = a {1 xt 5 (1 x) (a2)} 1, 
1 23) 
® 
ier: 4 R 1\f/1-v2V/ Fro (UHIe11) 
‘To I-v" TRY 3(Ay?—‘\4 E 2 
a (FE) + 3g 
2 eae _ 
3@ +) 
ene NE 2h wy? 
n l-v Oo, 5 R2 : Oo 


Ar = Ar ee sin 4w t 
O1- or Oo 


Wn 


Since the negative sequence magnetic field for a per-unit negative se- 


quence of I, = .1 is quite small, the maximum steady-state deflection Ar 
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should be negligible unless the ring natural frequency is very near 


"ww. 
oO 
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Chapter IV 
Electromecnanical and Thermal Effects of Faults 

The fault behavior of a superconducting generator as seen at its ter- 
minals (both electrical and mechanical) is essentially the same as for 
a conventional generator. The superconducting generator has the advantage 
of lower synchronous reactance, resulting in superior system stability. 
However, the effects internal to the machine produced by the large fault 
currents and forces may be more severe in the superconducting generator 
for the following reasons: 

(1) Alternating currents induced in the field winding could drive 
it normal. 

(2) Heating produced by induced currents in rotor parts can raise 
their temperature considerably, because of low thermal 
capacity at cryogenic temperatures. 

(3) Magnetic radial stresses upon the electrothermal shield can 
produce large deflections because of the lack of radial sup- 
port between shield and inner rotor structure. 

(4) The thermal isolations of the rotor make it structurally more 
vulnerable to large transient torques. 

None of these effects represents a fatal flaw of the superconducting 
generator concept. All can be overcome with proper design consideration 
of the effects. However, because of the critical competition between 
thermal isolation to minimize heat leak and structural integrity during 
faults, it is necessary to predict as well as possible the rotor stresses 
and heat loads during faults. Expressions for these stresses and heat 


loads are derived in Appendix III. In this chapter, the mechanical and 
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thermal model of the shell will be used to estimate the resulting motion 
of the shield, the temperature rise of the shield, and other determinantal 
effects upon the field winding. 
1) The Effects of Fault Torques 
The existence of large oscillating torques within a generator fol- 

lowing a fault has been known for many years and studied in numerous papers. 
The most complete analysis including initial loading of the generator for 
tnree-phase short, line-to-line short, and synchronizing out of phase is 
included in Ref. [7], Woodson's Report to Consolidated Edison. Appendix 
I presents a derivation showing the applicability of the results to 
superconducting generators taking into account the constraint on the 
various mutual inductances imposed by the presence of the conducting elec- 
trothermal shield. The cases considered are appropriate for the super- 
conducting machine, because the present armature configuration is a delta 
connection which precludes the possibility of line-neutral machine faults. 
Even when a machine is Y- connected, a large impedance is placed between 
neutral and ground to limit fault currents; the line-neutral fault can 
thereby be made less severe than the other faults in which such a simple 
current limitation is not possible. 

The calculation of fault torques for conventional machines usually 
results in a total torque acting on the stator (or rotor). Since both 
the field windings and damper bars or windings are set in the rotor iron, 
the distribution of actual forces between field windings, damper bars, 
and rotor iron is not so important. However, in the superconducting gen- 
erator, the damper winding is the electrjcal shield, which is a thin- 


walled tube supported only at the ends. The torque exerted upon this 
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tube can produce the following failure modes: shear off the shield if 
the yield seEees of the material is exceeded; cause torsional buckling of 
the shield; produce large torsional oscillations of the shield if the 
natural torsional frequency is near the driving frequency of the torques. 
Oscillating torques are exerted upon the field winding also, to a lesser 
extent. The field winding must be mounted with sufficient rigidity to 
prevent excessive motion of the superconducting wire. There has been 
some concern as to the possibility of driving the superconductor normal 
because of motion during faults. The analysis shows, however, that most 
of the torque is taken by the shield. This fact lessens concern about 
the superconductor, but makes prediction of shield stresses more critical. 
2) Torque Upon the Field Winding 

The field winding is subjected to an oscillating torque, as 

given by Eq. (A-42): 
- x! mi u 
r, - (3, *\ l a ae re (1, Yo, (IV-1) 
d d “d ety ac 
do 

Inspection of Eqs. [A-32(a) thru (e)] shows that var equals zero for 
faults from open circuit. Therefore, the alternating torque on the field 
winding is essentially zero for all shorts from open circuit. The second 
point to be observed from Eq. (IV-1) is that the factor [1 + (To 0 ou 4] 
represents an attenuation factor due to the presence of the damper wind- 
ing. For an electrical shield whose thickness is greater than a 
skin depth, the attenuation factor given by Eq. (A-53) must be used. 
The result is that, even for a loaded machine where an is not zero, the 


torque upon the field winding is considerably less than the total rotor 
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torque. Another way to state this is to say that a shield which is effec- 
tive in intercepting the asynchronous magnetic field seen on the rotor 
also intercepts the alternating torque which is exerted on the rotor. 
Equation (IV-1) is evaluated for the cases of three-phase fault, 
line-line fault, and synchronizing out of phase, and the results are 
listed on the following pages. Inspection of these results indicates 
that the per-unit torque on the field winding is not large for reasonable 
values of attenuation factor. 


Alternating Fault Torque Exerted Upon the Field Winding 


1) Three~Phase Fault from Load 
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2) Line-to-Line Fault from Load 
The line-line fault produces two oscillating torque components, 
one at line frequency and one at double line frequency. The attenuation 


factor must be evaluated separately for each. 
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Resistance-Limited Shield (IvV-4) 


x - x'\/ x" ia 
—_ (23) 2) (2a) aac 


b) Torque Upon the Electrical Shield 


For an effective shield, the field winding torque has been shown 
to be small. For calculation of the effects of torques on the elec- 
trical shield, it is then slightly conservative to assume that the 
total resultant oscillating torques, as given for various faults in Eqs. 
(A-33) through (A-37), are applied to the shield. 

If the shield is supported in torsion at both ends, the maximum 
torque along the shield,tending to shear it off azimuthally, is one-half 
the applied torque. If the natural frequency of the shell in torsion as 
given by Eq. (IV-3) is much higher than the oscillating torque frequency, 
the shield will yield in shear if the yield stress in shear of the mater- 
ial is exceeded. The frequency of the oscillating torques from Eq. (A-1-b) 
are fundamental line and second harmonic frequencies. The lowest tor- 


sional natural frequency is 


fo ele (IV-5) 
The quantity /G/P is roughly constant for most metals; evaluation of 
the equation shows that a shield nearly 10 meters long would be required 
to have a natural torsional frequency as low as 120 Hz. Even for the 
example 1000 MVA machine in section II-4, the shield length is only 
about 4.2 meters. The torsional natural frequency will almost always 


be well above the driving frequency of the oscillating torques. The 
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resultant state of stress in the shield must then be found by combining 
the shear stress due to torques with the centrifugal hoop stresses due 
to rotation of the shield. 

The maximum shear stress due to torque of magnitude tT for a shield 
of inner and outer radii R; and Ro is given by the following equation: 


Za T 
Oo 


(IV-6a) 
m(R* - R,’) 


The tensile stress due to rotation of the shield is given by the follow- 


ing equation: 


Pm WOR. 
2° = Dae (IV-6b) 


where oS is the mass density of the shield material. These stresses 
must be combined using Mohr's circle. For the ductile materials of 
which the shield will be composed, the most logical failure criterion 
is probably the maximum shear theory. The maximum shear stress result- 


ing from Oo. and 0,, above must be less than the yield stress in shear. 


iT 

The yield stress in shear is about one-half the yield stress in tension. 
The maximum torque Tt to be used in Eq. (IV-6a) depends on the fault 

condition considered. Expressions for these torques are given by Eqs. 

(A-33) through (A-37). The most severe which is of practical interest 

is the line-to-line fault at the terminals, while operating loaded 

at rated power factor. The case of over-excited operation at zero power 


factor at rated reactive power is slightly worse, but is a more unusual 


condition of loading. 
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A second possible mode of failure is buckling of the shell due 
to torsion. The situation in which this mode of failure can occur before 
simple failure in shear is that for a very thin cylinder. Most calcula- 
tions have shown that the critical torque for buckling is well above that 
for shearing off the shield; however, the criterion is reproduced here 


from Ref. [9] for completeness. 


Buckling of a Shell in Torsion 
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Unidirectional Shield Torques 


Unidirectional torques or induction motor type torques are exerted 
upon the electrical shield during transients. In general, they are 
small compared to the peak values of the alternating torques. They 
can be evaluated from the power dissipation in the shield, due to the 
asynchronously rotating magnetic fields. The heating in the shield is 
itself a more important effect than the unidirectional torques, and is 
treated in section IV-3. The unidirectional torque can be evaluated 


from the value of P ac? the power dissipated due to alternating currents 
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induced in the shield. Values of Ee for the various fault conditions 


are derived in IV-3. 


Torque unidirectional = = P 
O 


ac 





82 
IV-2 The Effects of Normal Mechanical Stresses of Electrical Origin 

In Appendix I it is pointed out that, during an electrical transient, 
when a magnetic field is established external to a conducting shield in a 
time short compared to the shield time constant, a normal stress is exerted 
upon the shield. Inside the superconducting generator, this field is dis- 
tributed sinusoidally in azimuth, producing a normal stress distribution 
which tends to flatten the shield. The magnitude of this stress is deter- 
mined primarily by the machine reactances as shown in Appendix I. However, 
the duration of this transient is dependent on the time constant of cur- 
rent decay in the shield, and the maximum deflection of the shield will 
depend upon its ring natural frequencies as derived in Appendix I. 

In this section will be presented expressions for the maximum radial 
deflection of the electrical shield under various fault conditions. 

Three~Phase Fault 

The normal stress resulting from a three-pnase fault from open cir- 
cuit is given by Eq. (A-61). In order to determine the flattening of the 
shield, the transient radial motion of a point on it can be solved as for 
a lumped-spring mass system subjected to a driving force whose time depen- 
dence is given by the normal stress equation. Equation (A-66) gives the 
radial deflection of a shell subjected to normal stress given by Eq. (A-65). 
The natural frequency of the shell is given by Eq. (A-65); with these ex- 
pressions an equivalent spring constant and mass can be defined for 
radial motion of the shell at angle y = 0. Evaluation of the second spa- 


tial harmonic terms of Eq. (A-61) gives the following stress distribution: 
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If we consider the deflection at y = 0 as a function of time, the terms 
containing sin 2y will produce no contribution; the cos 2y will produce 
a deflection as given by Eq. (A-66). To solve for the transient shell 


deflection, the solution to the equation is given by 


mx + kx = £(t') 


where f(t') is the time dependence of the cos 2y terms in Eq. (1V-8). 


Stress producing deflection at y = 0 is expressed as 
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From Eqs. (A-65) and (A-66), we obtain the following: 
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The complete solution of the differential equation with the driving func- 
tion (V-9) requires much algebra. However, the result being sought is the 
maximum deflection of the shield during the transient. To obtain that re- 
sult, we can make use of the fact that both the armature time constant and 
the subtransient time constants are at least several cycles of line fre- 
quency, and for most machines, even longer. Using this assumption, the 


differential equation can be written 


Ee fe 5 6 
mx + Kx = aiid ZF «cos 2u,t - 57 «cos wot 
fo Als 


2 io 2: 
The solution for x___ depends on the relationship of the shell's natural 
frequency oO = / K/m , aS given by Eq. (A-67), and the driving frequen- 
cies W and Zu) Damping of these radial motions has been neglected; 
therefore, the deflection for W = Wor w = 2W will become infinite. 
However, the limits for Ww < W and Ww > 2u can be estimated by the 


undamped case. The results are as follows: 
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From Eq. (1V-11), it appears that the deflection is less for a shield 
with natural frequency below W5 however, this corresponds to a shield 
with a low value of K producing a larger Xo° The stiffer high frequency 
shield produces dite smaller deflection. 
3) Rotor Heating During Faults 

The effects of rotor heating during faults was studied extensively 
in the past [22] and criteria of 1 *t for a generator was developed. The 
pt is essentially a measure of the amount of rotor heating due to trans- 
ient negative sequence currents which can be tolerated without permanent 
damage to the machine. In a superconducting machine, this rotor heating 


effect must be accurately predicted for the following reasons: 


1) Small specific heat of metals at cryogenic temperatures; 

2) Relatively small thermal mass subjected to heating; 

3) Possibility of driving the superconductor normal. 
One of the most important functions of the electrical shield is to inter- 
cept the large negative sequence magnetic fields during faults. Since 
it is a relatively small thermal mass, considerable temperature rise can 
be experienced in the shield. Attempts to add thermal mass to the shield 
are difficult; materials such as stainless steel with low electrical con- 
ductivity also have low thermal conductivity; during the period of the 
transient, very little heat is conducted into them. A thicker shield d 
high electrical and thermal conductivity adds effective thermal mass, but 
increases the shield time constant and reduces the shield's effectiveness 
as a damper winding, as shown in Ref. [5]. 

The calculation of fault heating employs the results of Appendix I 


to obtain the armature fault currents, the magnetic field model to obtain 
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the heating, and the thermal model to calculate temperature rise of the 
shield, due to the nonlinear temperature dependence of the shield's 
electrical conductivity and specific heat. 

Heating produced in the shield during faults results from the fol- 
lowing sources: 
1) The de currents induced in the shell to maintain its flux 
Linkage constant during the subtransient period; 
2) Alternating currents induced by the trapped flux linkages 
in the armature which appear to rotate backward relative to 
the rotor at angular frequency Wo} 
3) Alternating currents induced by negative sequence fault cur- 
rents resulting from unsymmetrical faults. 
The most severe case is the line-line fault, which results in wm balanced 
currents even after the transient period. The results for this case can 
be used to determine the temperature rise occurring before the system cir- 
cuit breakers interrupt the fault for a short-term fault, or the ultimate 
temperature rise for a sustained fault. 
Three-Phase Fault Shield Heating 
The magnitude of the increase in positive sequence armature cur- 
rent is given in Table IV-l1. The resulting current induced in the shield 
to maintain constant flux linkage at the shield is derived from the mag- 
netic field solution given in Appendix I. The resulting expression for 


power dissipation due to de shell currents is 
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qi? a from Table IV-l. 


The three~phase fault is symmetrical and no negative sequence armature 
currents are produced. The ac power loss in the shield results entirely 


from the trapped armature flux. 
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The total dissipated power is the sum of De and Pact 
In calculating the shield temperature rise due to the fault, a 
step~by-step calculation is performed,taking account of the integrated 
power dissipated and the dependence of conductivity and specific heat on 
temperature. It can be shown that the power dissipated in a three-phase 


fault is independent of initial current; i.e., it depends only upon 


terminal voltage before the fault. 
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Line-to-Line Fault Shield Heating 
The de power for the line-line fault is given in Eq. (IV-13) if 
the appropriate time constants are used. 
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The ac power is computed using Eq. (IV-14), using the appropriate 


value of magnetic flux density and aneular frequency. 


1 Gm R 2 
bee < 
277 R? 1 (wT) 
ts ee ee = 


Equation (IV-16) must be evaluated for both trapped armature flux and neg- 


ative sequence armature currents: 
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Trapped Flux: 


Negative Sequence (2 components): 


1 Xe 
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IV-4 The Field Current During Faults 
During the period immediately following a fault, the field current 


level rises due to the demagnetizing effect of the armature currents, and 
an alternating component of current is induced due to trapped armature 
flux and unbalanced conditions of the fault. 

The rise in the direct current level is not affected by the presence 
of the electrical shield; the rate of rise will, however, depend on the 
electrical shield time constant. The peak direct current level depends 
only on the increase in the direct-axis current during the transient 
period, that is, when the shield currents have decayed, but the field 


winding is still linking a constant flux. In per-unit notation, 


Ai = (x - x') Ai,' 
fae da} 
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The expressions for Ala are summarized for three-phase and line-to- 
line fault in Table IV-1. The largest rise occurs for a line-to-line fault 
on a generator operating at zero power factor over-excited. 

The alternating current induced in the field winding depends upon 
the attenuation of asynchronous magnetic field provided by the electrical 
shield. The equations are given as Eqs. (A-19) relating alternating field 
current and alternating direct-axis current. They are repeated below in 


per-unit notation, giving peak values. 
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The expressions for i, are given in Table IV-2 for three-phase short 
ac 
and line-to-line short circuit where i, is given by i, for a three- 
ac 


phase short circuit, and i, and i, for a line-to-line short circuit. 


Table IV-1 


Increase in Positive Sequence Direct-Axis Armature Current 


Following Fault 


1) Three-Phase Short Circuit(same for loaded or unloaded) 


V x" Cc x!" x'! - Tt xc 
* = " Sa a4 2. 
eer, apo ea re te 
“e 


dj x d 





os 


ve 
a = at ? 
A ay +r T) << << Ta 
2) Line-Line Short Circuit: 
(a) From Open Circuit: 
t-t, t-¢ 
ne 2) ee ow \ 
di 2x" x! + x i dk " (2 x! x" oa -| i dg 
a 2 
x" + x 
ae 
aa! oe 1B] 
Tah 2x" mer 
? Ve 
° as te t 
Mar” Xe xt TaeS< © << Tay 


(b) From Load: 


ve uf 
<r 1 ( A. Mie oua'e E. F 
di 2(x +x.) x" x + x. f “d2) 


x - x" 
a end Veeco’ Fat ; Veamcog s,} Be are 


Fue,’ F 320 given in Eq. (III-24). 


Au, © il + x, ! x" +x. - 
1 (x * x" sre ee. ooe 5.) 


- I, sin(6 + Op) 


09 
for t << Tae 





92 
- °. i 2x" + x. x + x, 2x"(x' + x.) : 
di 2(x"+ x.) x x + x, x" (x'+ x.)+ x" (x" + x) f 


eae 2x"(x" + x.) ae (x+x_) (x! - x") ) 


" (3 + x. low +x )+ x! carat ee 2 x_) (x - x") 


V cos 6, -~ V cor 6) — 1 sin(6 + 6p) 
sa S sa S a 


tt ’ 
for Ta << t << Tae 





1 2x'"'+ x x"+x 2x''(x + x_) 
Eg (22 = ( : ey, 
dj 2(x'+x_) x x + X, x" (x+x_)+ x(x" + x.) f£ Sa 


S 


cos 5} - 1, sin(6 + 6p) 


t 
for t >> Tag 





1) 


2) 


| 
| 


93 
Magnitude of Direct-Axis Alternating Current Following Fault 


Three-Phase Short Circuit (same for loaded or unloaded): 
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CHAPTER V 
Experiments 


Several experiments were performed during the course of this study 
to confirm intermediate results of the analysis. Some experiments were 
also performed upon the existing superconducting generator. These were to 
confirm prediction of performance, and to observe whether or not any unex- 
pected developments should occur. 

The first experiments were done to confirm predictions of armature- 
produced flux density, time constants of conducting shields, attenuation 
factors of shields, and power dissipation in shields. The experiments on 
the existing machine were to record performance for a sudden three- 
phase and line-to-line short circuit. Also, the generator was run as an 
induction motor with the shield acting as the rotor winding, to measure 
shield dissipation and to observe the effect on the field winding. The 
results of these experiments with explanations of the setup and comments 
on the results are given in the sections of this chapter. 

V-1 Armature-Produced Magnetic Field 

The magnetic flux density produced by the armature of the 80 KVA gen- 
erator was measured while the rotor had been removed for repairs. The 
armature currents were produced by a bank of three-phase transformers con- 
nected through an autotransformer to the line. The rms flux density was 
measured on the axis of the machine in the active region, and through the 
end-turn region. The results are plotted in Fig. V-l. 

From Eq. (A-48), the predicted rms flux density using the distri- 
buted armature current model was 140.5 gauss for armature current of 200 


amps. Using the discrete bar model with the Fourier analysis from Eqs. 
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(III-5) and (III-6) plus Eq. (4-48), the predicted value was 141 gauss. 


The measured value at the center of the active region was 139 gauss. 
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V-2 Shield Time Constants 

A number of experiments were performed to confirm the expressions 
for shield time constants which are developed in section A-2-C. Conduct- 
ing shields were used, whose thicknesses were about 30% of their outer 
radius, so that the diffusion time constants, as well as the lowest shield 
time constant, could be observed. The shield conductivities were first 
measured with dc, and tne time constants predicted from the conductivity 
and dimensions, using Eqs. (A-60,60'). The experiments were performed by 
placing the shield inside a motor stator and suspending a flux probe 
inside. 

Resistance was placed in series with the motor armature, to make the 
armature time constant short compared to the shield time constants to be 
measured. A step of voltage was applied to the armature and the result- 
ing voltage upon the flux probe was recorded. The recorded probe voltage 
was plotted on semilog paper, to determine the time constants. The 


shield time constant Ts and the first diffusion time constant T 52 could 


1 
be determined from the plots. The predicted and measured results are 
given in Table V-1. The measured values are all slightly smaller than the 
predicted values for the thicker shields. There are three reasons for 
this: first, the shield was slightly shorter than the overall length of 
the motor armature, therefore the two-dimensional analysis does not 

take into account the resistance contributed by the path of the currents 
around the ends of the shield. Second, the iron outer boundary of the 
motor had teeth for the armature windings; the inner radius of the teeth 


was used for the iron boundary in the equations, but the effective bound- 


ary is slightly larger, and would predict a slightly lower time constant. 





Third, it was difficult to make the armature time constant too much less 
than Too because the resulting flux inside the shield became very small. 
Although the armature time constant is too small to separate by the plot- 
ting technique, it provides an additional delay in the rise of the flux 


inside the shell, and appears to make observed time constants slightly 


larger. 
Table V-l 
Ts, (ms) Ts, (ms) 
Shield Ro(in) Ain) Oo Pred. Meas. Pred. Meas. 
Al 1.75 .5 2.28x10’ 7.0 6.8 47 58 
A1(N 1.75 .5 5.67x 10’ 17.7. 17.4 1.19 1.45 
temp) 2 
Bronze 1.75 a5 .65=x 10 2.0 1.92 ~134 ~139 
7 
Al 2.25 -116 2.11% 10 3.3 3.0 not measured 


eee <4 
_ “7 
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Figure V-2 
Shield Time Constant Experimental Setup 
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V-3 Attenuation Factor of Thick and Thin Shields 

An aluminum shield with dimensions 3.5 inches OD and 2.5 inches ID 
was placed inside a motor stator while the stator windings were excited 
with variable-frequency currents. The flux density inside the stator was 
first recorded as a function of frequency and current without the shield, 
The shield was placed inside the stator and the flux density measured in- 
side the stield at the same frequencies and currents as without the shield. 
The range of frequencies was such that the magnetic skin depth ran from 
much greater to much less than the shield thickness. 

The attenuation factor reciprocals predicted by Eq. (A-54) for a 
thin shield, Eq. (A-53) for a thick shield, Eq. (A-52') using the Bessel 
functions, and the measured values, are recorded in Table V-2. For skin 
depths greater than the physical thickness, the thin shield attentuaion 
factor is fairly accurate. For thickness greater than about two skin 


depths, the approximate expression given by Eq. (A-53) is quite accurate. 


Table V-2 
Ryeae 1.25" =e 31675 om, R= 44.45 mm, Y RR, = 37.6 mm, A = 12.7 mm 
WK ect 
= eo 
Bessel Thick Shell Thin Shell Exper. 
Function Approx. Approx. 
f 6 (mm) Eq. A-52' Eq. A-53 Eq. A-54 
40 16.67 2.46 leak 1.33 1.75 
80 11.78 4.52 4.68 3.56 3.8 
100 10.54 5.65 Dias 2 4.42 4.5 
200 7.45 12.20 12.14 8.66 oad 
500 4.72 47.2 47.6 21.42 47.0 


1000 30 196.7 197.6 43.0 200.0 





Figure V-2A 
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V-4 Power Dissipation in a Shield 


The power dissipated in a conducting, stationary shield inside a 
three-phase induction motor stator excited by balanced, three-phase cur- 
rents was measured. The flux density as a function of stator current was 
measured without the shield and armature i*r losses were recorded. The 
shield was placed inside the stator and the increase in power at the 
Stator terminals recorded for the same armature current as without the 
shield. The time constant of the shield had been determined independently 
by the experiment of section V-2, but was also checked by recording the 
attenuation factor of flux inside the shield. 

The flux density without the shield was recorded as a function of 
axial position through the end-turn region in order to determine the ef- 
fective length of the stator. The power dissipation was predicted from 


Eq. (A-50), using the flux density without the shield, the shield time 
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constant, and the geometric factors to be 191 watts. The measured dissi- 
pation was 196 watts. 
V-5 Sudden Short Circuits of MIT Experimental Generator 

The MIT first-generation superconducting machine (nominal rating, 
86.5 KVA) was subjected to sudden three-phase and line-to-line short cir- 
cuits from open circuit. Oscillograms of the resulting currents were 
recorded in order to check the expressions for transient and subtransient 
reactance, and for subtransient time constant. 

Several runs were made of the three-phase short circuit test. All 
gave about the same value for reactances, but the later values of sub- 
transient time constant were longer. This, probably, is because the shield 
was below nitrogen temperature for the later runs, since nitrogen was not 
used after the initial cool down. Tests made in September, 1969, were 
done at about 3% of what was then considered rated field current; results 
of the three-phase tests are stated below. The impedance bases results 
from the nominal machine rating of 64 volts, line to neutral, phase cur- 
rent of 450 amps, and power rating of 86.5 KVA. The theoretical impedances 
are predicted from reference [6], and the time constants for equations in 


Appendix If. 


Table V-3 
Quantity Theory Measured (per unit) 
x -0133 .0865 
x' .0188 2 .0104 2 -0675 
x" -0054 2 .0070 Q .0455 
Ty -0116 sec. 015 sec. 
Ai, -061 amp. -063 amp, 
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FIGURE V-3 
Oscillograms of Sudden Short Circuits upon 80 KVA Generator 
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The field current rise was measured for the three-phase short cir- 
cuit and was found to be popetetent with the measured reactances and the 
equation for field current surge in Chapter II. 

Later experiments have shown that a more realistic voltage rating 
for the machine is 45 volts, line to neutral, at a field current of 56 
amps. Another set of sudden short circuit tests was run with the field 
current at 6 amps or 10.7% of rated voltage. These tests were run with 
similar results. No adverse effect upon the machine was observed, and 


Ssuperconductivity was maintained. A line-to-line short circuit from the 


same open circuit voltage was also performed, with no adverse effects. 


V-6 Induction Motor Loss Experiment 


The 80 KVA superconducting generator was run as an induction motor 
by driving the armature through a three-phase autotransformer and a set 
of specially built step-down transformers. The electrothermal shield 
served as the induction motor rotor. The power into the armature termi- 
nals was measured. From knowledge of the armature losses and the bearing 
and windage losses, the power dissipation in the shield can be calculated. 
The shield time constant was measured before the run was started, and 
after the last point was taken, by imposing a voltage step on the armature 
and recording the open-circuit field windinp voltage. The shield time 
constant had changed from beginning of the run to the end, because of 
the dissipation, but at least the conditions for the first and last point 
of the run were determined. From knowledge of the shield time constant, 
armature currents, and rotor speed, the power dissipation in the shield 
could be predicted from Eq. (A-50). The experimental and predicted re- 


sults are given below. 
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RPM Ib I _ (amp) P roemeas « evesees Swiield ‘eiigeld (A-50) 
2400 O12 160 53 ZZ OL B21 
600 0062 160 115 15 100 90.0 


The machine was run to 3,550 rpm as an induction motor, voltage was 
applied to the field winding, and the machine was synchronized with the 
60 Hz. system. Armature current as a function of field current was then 
recorded with the machine running as a synchronous condensor. The result- 


ing V curve of armature current vs field current is shown below. 
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FIGURE V-4 


"vy" Curve for 80 KVA Generator 


8 Z 16 
Field Current (amps) 
Terminal Voltage = 9.7 volts line-line 


Rated Armature Current = 450 amps 
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V-7 Slotted Electrical Shield 

The shield time constant of a 
conducting shield can be decreased by 
cutting circumferential slots as shown 
in the figure at the right. The resis- 
tance to axial current flow is increased 
without affecting the inductance, thus 


the L/R time constant is decreased. 


UCL 


This technique permits a shield designer 

to use a thick shield for thermal capacity 

while keeping the time constant to a low Slotted Electrical 
enough value to achieve optimum damping, Shield 

as peplained in Chapter II. 


The cross section of one side of the shield is shown below, giving 


definitions of the parametersof the slots as follows: 


w = shield thickness 
s = slot spacing 


h = slot width 





d = slot depth 
t = w-d 
The effective thickness A, of the shield for use in calculating the 


shield time constant is given by 


e S 
h/w + 2 dn(1+5> 


The time constant of the slotted shield can then be calculated using Eq. 


(A-56) with A. substituted for A. 
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Experiments were performed using conducting paper with different 


numbers of slots and spacing. The resistance increase in the paper from 


one end to the other was measured as the slots were made. 


the experiment are given below, comparing the resistance predicted by 


the formulae and the measured values. 


Predicted Resistance 


Number of slots KQ/inch 
1 8.69 
Z 10.80 
3 12.60 
5 14.55 
9 16.33 


17 17.68 


Measured Resistance 


KQ/inch 


8.65 
10.60 
12.30 
14.70 
16.50 
17.90 


The results of 
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CHAPTER VI 
Conclusions and Suggestions for Further Study 


The purpose of this study is to identify the critical considera- 
tions in the design of a superconducting generator rotor, to insure that 
it can survive fault conditions. The classes of fault conditions and 
criteria for survival are as follows: 

1) A short-term fault of from five to 15 cycles which is 

cleared by the system circuit breaker. 
Criterion: Superconductivity must be maintained. 

2) A sustained fault inside the system circuit breaker which 

cannot be cleared. 
Criterion: Structural integrity of the rotor must be 
maintained such that permanent damage does 
not result. 
Chapter II has presented a discussion of the critical effects with a pro- 
posed design for a 1000 MVA machine which satisfies the above criteria, 
while maintaining satisfactory steady-state performance with the inherent 
advantages of a superconducting generator. Chapter III, IV and V have 
presented the analytical tools and experimental verification used to 
develop the design in a form which can be used to evaluate fault effects 
upon future designs. The conclusions which can be drawn from this study 
concerning the critical considerations in design of a superconducting 
generator to survive fault conditions are presented in the following 


section. 
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VI-1 Conclusions 


1) 


To insure the maintenance of superconductivity through a short- 


term fault, the following design considerations will be critical: 


(a) 


(b) 


(c) 


2) 


An attenuation factor no greater than .01 at line epequeney 
must be provided to limit temperature rise in the 4.2° K 
region. 

The superconductor operating point must be established with 
anticipation of the field current rise which results during 
the period prior to breaker interruption of the fault. 

It is possible to provide a shield which will produce nearly 
optimum damping consistent with the attenuation requirement 
of (a). The dynamic stability should therefore be superior 
to that of a conventional machine, where it is more difficult 


to tailor the damping coefficient. 


To maintain structural integrity of the rotor, the following 


design considerations are established: 


(a) 


(b) 


For a shield with an adequate attenuation factor, most of 

the transient torques are taken by the shield. To provide 
adequate support for the shield, it should be run at room 

temperature. A cryogenic shield with adequate structural 

support will have an intolerable heat leak. 

The thermal capacity of the shield can be made large com- 

pared to the thermal capacity of conventional machines of 

comparable rating. Values of te up to 20 seconds are 


possible without excessive shield temperature rise. 





Pi 


(c) Provision for shield structure to support torques will also 
provide sufficient support for transient normal mechanical 
stresses upon the shield. 

(d) Mechanical natural frequencies of the rotor will be high 
enough not to cause problems for a shield with sufficient 
structure to withstand the transient mechanical stresses. 

In summary, it is possible to design a rotor for a 1000 MVA super- 
conducting generator which can survive short-term faults without loss of 
Ssuperconductivity and sustain faults without damage. This can be done 
while retaining the advantages of lower weight, lower synchronous reac- 
tance, improved efficiency, better damping characteristics, larger Isat 
and higher rotor natural frequencies than for a comparable conventional 
machine. 

VI~2 Suggestions for Further Study 

1) The most important question which is yet to be answered is 
exactly how to relate short sample characteristics of supercon- 
ducting wire to the transient conditions in the field winding. 
Just how much alternating current, what rate of current rise, and 
how much of a temperature increase will drive the winding normal, 
must be determined experimentally. With this kind of information, 
tne analytical tools developed in this study can be used to arrive 
at a shield design which will permit best use of the superconduct- 
ing winding. 

2) A restudy of weight optimization should be done, taking into 

account the rotor structure necessary to support the transient 


electromechanical stresses. 





EZ 


3) A complete study of the refrigeration and shield cooling re- 
quirements should be done. It may turn out that the optimum shield 
temperature is somewhat below room temperature, if the rotor exit 
gas is used to cool the shield. Otherwise, conventional gas cool- 
ing of the shield at room temperature would probably be best. 

4) The specifications for the secondary or thermal shield should 
be more exactly determined, once its optimum operating tempera- 


ture is determined. 





113 


A-1 The Lumped-Parameter Machine Model 


Generators are conventionally mod- 





eled as shown in Fig. A-l, where each 
winding has a lumped inductance and 
resistance plus mutual inductance to 
the other windings. Three identi- ‘ 
cal windings 120° apart in space rep- 
resent the stationary armature, two 
orthogonal moving windings represent 
the damper windings, and a moving 
winding represents the field winding. 


Use of such a model whith appropriate 


engineering approximations and mathe- 


Figure A-l 


matical transformations yields 
expressions for transient reactances and time constants describing the 
behavior of a machine following a fault. 

A superconducting generator can likewise be modeled as in Fig. A-l. 
However, the damper windings for the superconducting machine consist of the 
electrothermal shield. During transients, the process of magnetic diffu- 
Sion through the E.T. shield (and the stainless steel dewar shells and 
metal structure), as well as the decay of currents in the shield's equi- 
valent IrRcircuit, affects the observed behavior. The adequancy of rep- 
resenting the electrothermal shield as a lumped L-R circuit, i.e., as one 
time constant, is investigated in section A-2-c. The result of that anal- 
ysis shows that the first diffusion time constant is less than the decay 


time constant by a factor of approximately A/R, i.e., the thickness-to- 
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radius ratio of the E.T. shield. For the thickness-to-radius ratios 
being used, the determination of stator currents and torques following 
faults will be sufficiently modeled by neglecting the diffusion time 
constants. 

The ultimate aim of this investigation is to determine the electrome- 
chanical and thermal effects within the rotor of a superconducting gen- 
erator following a fault. The purpose for using the lumped-parameter 
model is to determine the magnitudes of armature currents following vari-~ 
ous types of faults based upon the values of reactances analogous to those 
for conventional machines, but defined appropriately for superconducting 
machines. Out of this same analysis comes magnitudes of total fault torques 
exerted upon the armature. The total rotor torques must, of course, be 
equal and opposite, but their distribution between electrothermal shield 
and field winding must also be determined. The resulting armature fault 
currents are used in the magnetic field model, as developed in A-2, to 
determine the electrical forces and electrical heat loads to which the 
electrothermal shield and the field winding are subjected during fault 
conditions. 

A-l1(a) The Lumped-Parameter Model Including Resistances 
A lumped-parameter model as shown in Fig. A-1l can be used in 
two ways. If the resistances are included in each winding, the differen- 
tial equations describing the circuit behavior, even after appropriate 
transformations, can be solved correctly only by using a step-by-step 
numerical integration with a digital computer. Such investigations are 
necessary for system performances and stability studies, and have been 


carried out in reference [5]. However, for determination of internal 
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effects within the machine itself, it is felt that more insight into the 
important aspects of fault behavior can be obtained by use of approximate 
analytical solutions. 

The approximate solution of a six-winding model of a generator for 
fault currents and torques is by no means unique. In fact, the Park's 
transformation idea which led to many papers on the subject was published 
over 40 years ago. The one unique aspect of applying such analyses to 
superconducting generators results from the absence of iron in the rotor; 
that is, the fact that the generator is an air-core machine. Convention- 
al analyses have normally assumed that all windings of a common axis are 
coupled by a common mutual flux. In an air-core machine, this is obvious- 
ly not true. However, the physical placement of the E.T. shield and its 
continuum nature lead to another simplifying assumption. This is that all 
flux wnich links the armature phases and the field winding must also link 
the E.T. shield (or the damper winding representing the E.T. shield for 
the lumped-parameter model). This assumption leads directly to a rela- 
tionship between the mutual and self-inductances of the model shown in 
Fig. A-l, and is stated algebraically in Eq. (A-10). 

The assumption produces a simplification similar to the common 
mutual flux assumption for conventional machines, and makes an approximate 
solution of the machine equations possible. This assumption was used by 
Einstein [5] in preparation of his computer model, and is demonstrated 
alpebraically by Kirtley [6] from a magnetic fields model to determine 
mutual and self-inductance. 

The effect of resistances in the lumped-parameter model is to cause a 


decay of fault currents; the mapnitudes of the fault currents are 
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practically independent of the resistances. Therefore, once the appro- 
priate time constants for decay of currents are determined, the model 
excluding resistances can be used. This provides considerable simplifi- 
cation, especially for non-symmetrical faults which are in some ways most 
severe upon the machine. 

In the present section, the resistances are included to determine the 
appropriate definitions of the transient reactances, time constants, and 
attenuation of asynchronous fluxes. 


The lumped-parameter 





model for the superconducting 
generator is shown in Fig. 
A-2. The E.T. shield is 
modeled as two R-L cir- 
cuits, one in the d 

axis and one in the q 
axis. A general rela- 
tionship between the 
fluxes linked by the 
windings and the cur- 


rents is given by Eq. (A-1). 


Figure A-2 





LT 


r ars ard Trl M M M i 


a aa ab ac af akd akq a 

rb Stemeaue ose “oe “tea “Sika it, 

oc : oe a Mob - mee coe Med M kq i. 

Ye | | Meg — Men ~ Meg Le Mexg 9 te (A-1) 
ara ie dh ede Tere kad tid 


ae fr Mega y aap 7 Mac y : Lig ce 


Most of the inductance coefficients are functions of the angular po- 
sition of the rotor 9. Use of the Park's transformation will yield a 
considerable simplification of the matrix. However, first the number of 
coefficients can be reduced by reciprocity and symmetry. Reciprocity 


implies that the matrix must be symmetric about the diagonal; for example, 


IM! = D,! 
Magl = [Mega 


etc. 


The stator self- and mutual inductances will be a constant because of no 
rotor saliency. By symmetry, the magnitudes of these inductances must be 
equal. The rotor-to-stator mutuals must be dependent on the rotor angular 
position. As in conventional machine analysis, only the first spatial 
harmonic of flux will be considered in the lumped-parameter representation. 
All angular dependences must be sinusoidal. By physical symmetry, the 
magnitudes of the sinusoidal dependences must be equal for similar sets 
of windings. Mutual and self-inductances within the rotor are constant. 


By symmetry 





118 


Stator 
a . a . oe i 
one - ae = pe = aa 


Stator-to-Rotor 


Moxal * Moxa =lMygl = ete. = My 


Mie] = [Mel = [Ml = Meg 
Rotor 

Maed 7 ek 

Bee ka > 


Park's transformation from a,b,c variables to dqo variables is defined 


below for the dummy variable, u: 
27 27 
2/3 cos 8 2/3 cos (6 - 3°) 2/3 cos(® + 3? u. 
27 27 
u ={!- 2/3 sin 9 =- 2/3 sin(@ -->) ~ 2/3 sin(® + => U) 


Als 1/3 1/3 1/3 u 


When applied to Eq. (A-l1), the result is as follows: 





fd d 
0 “La 0 0 0 Meg " 
0 0 J 0 0 0 i 
O O 
-3/2 Mea 0 0 be Me. 0 te 
2 i 
3/2 Mg 0 0 Mey L. 0 kd 
0 -3/2 Med 0 0 0 Li. kq 
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The inductances L Lib? Le» Le» Oren Meg? and M are defined in terms 


fk 
of the geometric parameters of the superconducting generator in Ref. [13]. 


The voltage equations associated with the dqo transformed flux equa- 


tions are well known, and are given below. 


direct axis 


ao rig + pA, = hac 
Magee o = tle TP Aga 
ee * rei- +p de 


quadrature axis 


Le = = Pan + PA, + WAG 
vie = QO = Beane +p She 
p = d/dt 


The next step in any machine analysis of this type is to convert to a 
dimensionless set of parameters called the "per unit system". Base 
values of current and voltage are chosen, usually the rated machine values 
in the armature, and the corresponding base values for the other windings 
are derived. In determining base currents for rotor circuits, account 
must be taken of the fact that the direct and quadrature axis armature cir- 
cuits are 3/2 times as effective in producing armature flux as the rotor 
circuits, even when the turns ratio has been accounted for. This is lecause 
a three-winding armature is now being represented by two windings in the 
d and q axes. The derivation is carried out in most machine theory texts 
and is also given in Ref. [5]. The exact expressions for the base values 


of flux and inductance are not important to the results here. However, 
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after the appropriate normalization has been done, the per unit flux and 


voltage equations can be written as follows: 


Per Unit Machine Equations 





flux 
Na we edie de aa 
‘ae | Max Ly Mee Lok (A-2) 
re eae Ee Le 
_ aes M i 
g d dk q (A-3) 
aes ke nk A ak 
voltage 
wh 
v, =-ri, +2 y.,-— j 
d ad W, d Wo q 
= P_ = 
aCe Tae w hak ° 
= eee) a 
ee oe ws | (A-4) 
Pp “m 
a S — roi + h +—) 
q a sot) 6 CU 
< Pi * 
“4 Tle + w Ava 0 


All voltage, current, and flux quantities in Eqs. (A-2,3,4) are per unit, 
W is the angular line frequency, and Ww is the machine mechanical angu- 
lar frequency. For example, a value of i = 1.0 would produce rated 
open-circuit armature voltage. 

At this point, it would be useful to derive the constraint upon the 
mutual and self-inductances imposed by the presence of the E.T. shield. 
This result was referred to as the "shell constraint" in Reg. [5]. It 


will be derived differently here, but will yield the same result. 
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In Fig. A-2, we have placed three windings 
on a common axis. We assume that all flux 
linking windings f£ and d must also link 
winding k. An alternating current iy 2 Thy y 7 a 
imposed in winding d and the flux linked by 


d, i.e, the alternating open-circuit voltage 


upon d is measured. Figure A-2 
— 
d I, 


The amplitudes of the alternating quantities are given by da and I At 


dad’ 
the same time, we measure the flux linked by k, with k open circuit and f 


open circuit: 


‘i = Maly 


The fraction of flux produced by d which is linked by k is 


Since all flux produced in d and linking f must also link k, we can 


express the ratio of Ae to Ay Similarly to the above equation, 


£ fk 


In other words, since no flux produced external to k can link f without 
linking k, the effect is just as if the flux in k were produced by a current 


in k. We can solve the above equation for re in terms of i 


fea ; 
f i d 





rae 


But, by definition, re = Maela 


M., M 
Therefore, M qa : 


df Ly 


The relationship can also be expressed as 


Seep tee Moen ee 0 (A-5) 


The expression for C appears in later derivations and can be set equal to 
zero by the above assumption. 

Machine transients are conventionally studied in terms of subtransient 
and transient reactances. These are the reactances of the armature dur- 
ing the period during which the damper windings and the field windings main- 
tain constant flux linkages after a short in the armature. To determine 
the definitions of these quantities and definitions for the decay time 
constants of currents during these periods, Eqs. (A-2,4) will be solved 
for a symmetrical three-phase short circuit upon the armature. 


Using (A-2,3,4) we can write the following matrix of equations: 


~ bg Mae Mok 14 ra 
Pe. je 
Mago Te theo Mee te f= | SE 
oO oO 
2. P_ ip + Pp 
oe we Me ur rT hy we dak 0 


We can solve for i, using Cramer's Rule: 
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2 
ber + Lt. Leb Me 
To eee +# as) heen 
o f k oO fo k e M Po r, M 





























r 
eo f df o*k df 
d L ~ rb 
d D fd D 
(A-6) 
or I, wer | oM 
fre “kr¢ dfék “akéf 
D=1tpl\ Gg r.L ~ wW L.r.r 
oki f o dk f 
2 ——= 
_ Lely Mee Ma Mabe Mag g) A an 
ek KE 
Spay Dea tot eT k 
We define the time constants as follows: 
L 
_— - 
1 oS f 
1+ — 
2 Ws k 
Ly Me 
T = [1 - ] 
3 Work Bloat 
Le Mat 
. Wore {1 - rere (A-7) 
2 
Mak 
T = “k (1 -—— ] 
6 W Fr L 
ok d M2 
oo Psi, _ _dk ] 
ae Ly Lely Lila 
7 WoT} Mie 
fd 


Equation (A-6) can be solved for dg and rewritten, using (A-7) and C = 0: 


2 
ee CE, 
d pela Cae et dad reGn 15) p+ Tap re f 
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Inspection of Eq. (A-7) shows that 


6000 7 5 
a? 3 1 


since T» T. are related to the field winding time constant, and Ths Es 


Ts and T, are all related to the E.T. shield time constant; 


? 2 
~.e Lt (I. +T,)p + es a » % (1 + Top) (1 + Tp) 


and i Cech) pi a % (1 + T pd + TP). 


(1+ T.p)(1+T Jp) M 


d (+Tp)@+Tp) “dd” (+T p)@+T p) ry “f (A-9) 


From (A-9), the open-circuit and short-circuit time constants are ob- 


tained and listed below. 











be 
aaG - eI “7S r 
of 
2 
7 = T «= ak. {1 = Mee } 
do 3 eel L Leg 
2 
A-10 
d 5 Wr Lyle 
2 2 
— a Me ¢ ia - wk 
T) k L Li. 
Ty 1," M 
otk (1 - df } 
Le¢ly 


The general operational inductance can be written as follows: 


a Cl +T Pp) (1+ cys p) 
DP) ete ot. 7. (A-11) 
d d (1+To op). +T, Pp) 
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For times short compared to T,, i.e., 1/T, <p < ©, the inductances 


reduce to Li , the subtransient inductance. 
Ta 7 
n = BG 2G a 
d 3 TG0 T30 


Mea, 
Tai, 


Le L. (= (A-12) 


d d 


For times T, <t< T. pe ea =< p < ne » the operational inductance 
1 


T; 


becomes the transient inductance, La ; 


wh 
d 
bee, = 
d d Tao 
M 
df 
La = LQ “ime (A-13) 


In per-unit notation, inductances and reactances are numerically 
equal because the machine is operating at rated frequency, that is, for 


W@ = W = 1, The more conventional notation is to use per- 


0” “per unit 


unit reactances rather than inductances. 


Per-Unit Reactances 





d d M2 
tt ~_ = 2s) a 
Xi = Xqy (1 ey (A-14) 
Me 
df 
x. "ce" ye (] = ) 


More useful forms for the time constant expressions are in terms of the 


reactances as follows: 
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do £ 
x) - x5 
ee = T Gua Symmetrical fault 
O Si: \x, - Xy 
«' time constants (A-15) 
d 
T, = T= 
d c Xa 
x. 
ee = a? d 
ai Tao x4 


where Te is the field time constant, and TS is the E.T. 
shield decay time constant. 
For non-symmetrical (line-line) fault, one phase of the machine remains 


- Transformation of 


connected to the system through a system reactance, x, 


flux equations for this case shows that the effective direct axis armature 

inductance is increased by one-half the system inductance. Therefore, the 

same definitions of time constants can be used if, in Eqs.(A-15), the value 
of x, is replaced by Xa + = The resulting time constants for line- 


line short circuit written in form similar to Eqs. (A-15}are as follows: 


X47 Xq 
Tag = te |+- Xs 


cae 
x4 = x ~s + x Line-Line 
? se : ee ena =~ 
Tae ee x! = x" x Fault ome 
d d Ss ' 
+ x 
d 
i. we 
enema | yet 
qz S} Xo 
aa 


For transient following synchronizing of the machine to a system through 
a reactance x., the effective value of xX, must be replaced by Xs + Xo 


with the resulting time constants: 





T27 


1 


Te = oT i ace 
dc& f Xa + x 


s 
: + 
x x x4 : 
a, a i — d wo are ts Synchronizing (A-17) 
vd” at arbitrary 
ae mutt 
a 2 a xq | phase angle 
Z cL we 1 x, +x 
q d Ss 


Field Current During Transient 


The determination of distribution of transient torques between field 
winding and electrothermal shield as developed in section A-l-b requires 
that field current during transient be determined. This obviously can be 
determined only from the model including winding resistance, because of 
the shielding effect of the damper windings. This shielding effect is 
a function of only the damper winding time constant in the lumped- 
parameter model. For an electrothermal shield thick compared to an 
appropriate skin depth, the shielding expression must be determined from 
the magnetic fields model in section A-2-b. The lumped-parameter deter- 
mination of field current will indicate how the thick shield attenuation 
should be used for shields for which the thin shield, or single time 
constant attenuation, is no longer adequate. 

To determine the induced alternating field currents, we write the 


direct axis voltage and flux equations. 


oo p_ z 
Va Cae we da 


o” [se 
fe 


a Pp 
2 Tt an + we ak 
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f ft ee 

ra a aie Mela aa 
hac] = {7 Max be bee tale 
Ms ee a ae 


We can reduce the above six equations to one equation in Les 14 er 
uSing the same time constant approximations as in deriving the direct- 


axis operational impedance: 


i. = ot  { t)+ (2) i 
f ' '" r t " r d 
(L+T, p)Q+T, p) f (L+T, p)(l+Ty. p)\ £ 


To find the induced ac field current due to an alternating i,» we let 


i, = - I, cos wt and er = constant, and solve for 1, with the 
t 
> . 
assumption that Tao W > 1 
M I eee Oy le 
a neetis 2), oe vet! te 
t t 2 tt 
ac Tate vl +(T 342 Y1+(T, w,) 


(A-18) 


The square root factor in the denominator represents the damper attenua- 
tion. Further simplification results for the assumptions of resistance- 


and reactance-limited E.T. shields (equations in per-unit). 
Resistance-Limited E.T. Shield 


Ite | = (x, = x) | i. (A-19a) 





Leo 


Ore XK CX. =X.) 
- SS lt, | (A-19b) 


i 
| fe (x, - x, att Bee 


I 


A-1(b) _The Lumped-Parameter Model Neglecting Resistance 


The magnitudes of transient torques and currents are ffected 
very little by the resistance of the windings. By solving for these mag- 
nitudes in a model neglecting resistances, and applying appropriate decre- 
ment factors using time constants developed in the last section, the 
electromechanical effects on the rotor can be determined. In order to 
include the possibility of an unsymmetrical fault, an a -§8 transformation 
will be used instead of the dq transformation. Thea-§8 stator windings 
are the two windings equivalent to the three armature windings fixed in 
space, rather than the d-q windings which rotate with the rotor. 

The per-unit d-q inductance matrix given by Eqs. (A-2,3) can be trans- 
formed directly to the a, B notation. Equations (A-2,3) are rewritten 


below within one matrix. 


LX, a ig 
di ie a 0 ar 4, 
ae 0 8b, 6M (lO i. 
red eercgeeee ely Z a 
kq : Gea 2 , Ly tg 


The a,8 to dq and the dq to a,8 transformations are as follows: 


|| 8 sin | | ‘ - 9 cos 7 — 
i f : as 
A sin 86 cos 8 ig : ig cos 8 -sin 8 7 
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The dq inductance matrix in a,8 notation is as follows: 


A - Ly 0 M, 4cos o M 4cos y M, ySin ” a 
BB ° Bea M..sin 0 in 9 s 0 i 
fd aa de” B 
re a -M, 4c0s 8 “M, sin 9 Le Oe 0 i, 
ea -M, {C08 x3 -M, sin 9 Meg L 0 ing 
ave Moin Q ~M, {C08 9 0 0 Ly ikq 


Assume a short on phase 8 at t = 0, Ago = Meg (t=O) and solve for tear 


r 
5 oe koe y kd, 


Kd La lea 4 


We define the quantity Ln as follows: 


ce 

qi Leg kao (A-21) 
fis oats 

kd | kd 


We solve in Eq. (A-20) for Malka + M..i, with Go we h,(t=0), r 


fd f kdo 


dig t=) - 


"kd kd -( mz a i UM aMen fo ~ Meabercaot Mea Cheater Hea ta! 





1 
Mod --( 
oe Me” Uy be} (Me gla com MegMerARdot Mea Meghan MeaMeg ty! 


Since the electrothermal shield completely surrounds the field wind- 


ing and is a continuum of conducting material which permits currents to 
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flow such that its total flux linked must remain constant, a sudden change 
in direct axis current will produce no suddent change in ip: Therefore, 
the coefficient of i, must be zero. This is equivalent to the assumption 


d 
expressed in Eq. (A-5). 


ee ka tk “= 8 


The assumption that C = 0 leads to a more simple result for the sum of 


direct-axis rotor currents: 


ec kar + te. eee ea 
eee yf ee eee i 
Siigec fdf | oo ue | kdo * i d 
7 i a Sakae 
We define the quantity Ve as follows: 
“kd (A-22) 


an x L, eis 
Use of the constraint C = O and the definitions of io and 1g lead to the 
solution for and dg : 


a 
d 
Medtka * Meate 7 Ya, * TD 7a 


ea LQ - it, 


ne (A-23) 
= Lt Jig + We sin 90 + von cos 9 


————) iv + Ve cos 9 - ae sin 9 


he i L,@ - 


Equations (A-23) are now in a form which can be used to calculate the 
effects of unbalanced (or balanced) transients from loaded (or unloaded) 


conditions. 
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The machine configuration is as shown below with the a and 8 phase con- 


nected to an infinite bus Vag vase through a system reactance Xo 








Figure A~3 


The steady state of the system before the fault is represented by the 


phasor diagram below: 


8 = Wt 
Ve = Meale cos wt 
ae ev sin (wt-6) 


v =-=-v_, sin(wt-d6s) 


sa sd 
i = I, sin(wt- 6- 6p) 
i, = I. cos(wt- 6- Op) 


eh Ds sin(é6 + Op) 


ie I, cos( 6+ 6p) Figure A-4 
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Equations (A-23) can also be written using (A-12,13) and reactance notation, 


as (in per-unit): 


A Sse xi + Va, cos 6 - ven sin 6 


a da 
(A-24) 
= - xi, + ~ 
dg Xyig Wa sin 0 Vo cos 8 
and in dq components as 
= 2 i it 
Ya Yea t tabs Fy ail deur 
(A-25) 
= +iL = - x i + 
Yq "sq q s d “q Yar 
The torque acting on the stator is given by Eq. (A-26): 
i eg = Mesa 
= Cidade ea ¢)tq) 7) Seatwee ta (A-26) 


From the equations for Lig ene eedaid ui aS ecere we obtain in the quadrature 
axis the following equations: 
2 
Med 
+ L i 


se at 


= ven + (CL, - Lyi, 
In the direct axis, similar results are obtained. 
Med 
Medkd * Meats = Yay 7D *e 
Mey? ey Ep aa 
The torque equation can then be reduced to the following simple result: 


POY ita Venta (A-27) 
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The following equations relate torque angles, terminal voltage, system 


voltage, and voltage behind synchronous reactance. 


I xX cos 8p 
Se 
a I x sin Op 


6 = tan ? 


~1 Xe I, cos 9 
s V-x I sin Op 
C Sa 


1 
V. = x1, sin Op +[V.. - (x 1,)“(1 ~ sin@p)*] 


E2 = @. + xI, sin Opyr + (xI cos Op)? 


hy = sa + ix, =~ i + Vay cos OPT an sin 6p 


Na einer ode 


a ge Jeet + 
B 38 x aig + Yay sin Op ue cos 6p 


g%*5 


Eliminating hy? ig » we obtain the following expressions for Nae and 58° 


r = (x" as x) - ae cos Wt - ven sin wt 


sa 
(A~28) 
58 = — (x'' + x.) <- Va, sin wt + ee cos wt 
We solve for the currents in (A-28) with the following results: 
i.e - grid + grit cos ot ~ sin wt 
beg ¥ (A~29) 


1277 3x T + gif sin wt + gre cos Wt 
s 
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Line-Line Short Circuit from Load 
At t = 0, the 8 winding is shorted, trapping flux Yeo - Equations 


(A-29) become: 








WV p p 
= ese, ed ees eae 
(A-30) 
p p p 
d 
Le = = = + =i sin wt oe cos Wt 
The currents can also be expressed as 1 and oh 
te 
w y YW, (2x"+ x_) 
pesos rece ats — sin ut + =o 
d x + x, x 2x (x +x.) 
S 
TG + x) Yar cos 2wt - Vee sin 2 wt) (A-31) 


p p wy (2x"+x_) 
sa Bo } Ss 
q x + x Sin Wt vom cos wt + Px (xi x) 


Ss 
“Ltt. lf 
+ xx rx) (bay sin 2wt + ee cos 2 wt) 


We substitute the above into the torque equation: 


T= Ysa wy sin wt + cos Wt) 
sony Xe dy qi 


“Bo vy cos wt - pp sin wt) 
x"! dy qi 


Ss 


————— 2 2 
2x "+ x) a) a cos 2wt + a = Van Ysin 2wt] (A-32) 


Equations for Yar? en as functions of t following various faults were 
derived in Ref. [7] by Woodson; these expressions are reproduced on the 


following pages. 
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Line-to-Line Short From Load 


x" +x ue 


= ~ ee Oo A ta 
a x + x. Ueto ds, 55 x + x. en Sue 5, aa (A- 32-1) 
x= x" 
ean a. x. Ve. sin 6 Foe 
8 q Ss ve 
: 2 A 2x (x' + x.) . t/To 
dQ, xGerre ) + x(x" + x.) 
tteot 1" iS ' 
, 2x (x' + x.) 2x" (x + x.) : t/Tay 
[MC flO CSN eee 
x(x + x.) +x'(x + x.) x (x+x,_) + x(x +x.) 
2x"(x + x_) 
ae > ore 
x (x + x.) + x(x +x.) 
ae ttt i 1" 
: ; 2x"(x" + x.) ea i (x +x.) (x x") : t/Tao 
= Mpa of Ct” OUT Sn 
dQo2 x (x + x, )tx" (x saci) xX 2(x'+x_) (x - x") 
aL ‘oon We re 
7 2x" (x" + x.) aut : Caxctx_ ) Cx x'') x''(x x.) = 
i i ee aT 
x" Gtx. ) + x (x!"+x.) x-x" 2(x"+x_) (x-x") x" Gckx, )ta Cx", ) dt 
x"(x + x_) 
+ ee _ 
x (x +x) + x(x +x.) 
oe t | tf 
x(x + x.) 7 op x(x + x.) 
MS ee ci |, Ce Gk ye 
q x (x +x.) + x(x + x.) x (x + x.) +x (x + x.) 
Vc p ea: i 65) 
Veo = N.. sin (wt, = 5) Ag 
ee ie t/t.» 
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Line-Line Short from Open Circuit 


In the preceding, let = © 


re *B0%8 
x 





ae cos Wt Tia sin 2wt 


nee = ‘e sin WE 
- t/t 
An = e a g 
VARAR deco Udit 
2x" ~ Eta 2x"! 
F ato (1 ~ x) +x" : ae S ry yf 
i = 0 
a 


¥g 


1p == =F + Sh sin wt 


Three-Phase Short Circuit from Load 


Equation (A-32) for torque becomes: 


(A-32-2) 
a ' 
axe C/T ay aa" 
x+ x" x+x" 


The three-phase short circuit is independent of time of application 


of the short; therefore we apply it when flux linking a is zero - 





when wt = 6 + 7/2. 
oO 
Let y = Q 
sa 
YB , Me 
x = 0 
Ss 
From Equations (26), 
Vay p 
i = aoe cos Wt - ace sin Wt 
O x x 
Va A w 
i,=- aF + $1 sin wt + <1 cos ut. 
3 x! x 


i.e., 


(A-32-3) 





From Eq. (A-28) , 


VA 


From Ref. [7], 


oe eae eerdy 


ua, = = (x - x')I 


a 
_ cL Seer) © 


_ : x"(x- 
oc Ts x'(x- 
i. etter 

B Oe e 
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Wt - sin wt 
cos van n ) 


- (x - x )I. sin (6 + oe 


a CoS (6 + na 
t-t t-t 
. a x" x —— x" 
T =~ er A 
d + (Sr = Je d + zs 
es — t? = a ? 
x’) x' (x- x") 


a 


Three-Phase Short Circuit from Open Circuit at wt = 1/2 


I = Q before 
a 
= tt 
in ¥a,/% cos 
7. 2 a, Yaa 
B x'' x" 
VA 
c a x" Vay 


fault (A-32-4) 


Wt 


sin wt 


Var = M..I. F,. where M,,I = V 


fd fo di 


fd fo a 
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Closing at Arbitrary Phase from Open Circuit 


We assume the generator is out of phase with system by angle O.: 


a Meateo sin wt, Vg = Meateo cos wt (A-32-5) 
Yo = Meats, cos wt , Ye = Meats, sin wt Sanne 
a? Ne sin(wt - 65). oe ve cos (wt - 6.) Senin 
Yoo = View cos(wt - 6.) : V58 = Nee sin(wt- 6.) 

1 i ee aL 

BE “fo 


Closure occurs at t = QO. 


1 
i, = x" + x {ya = Vege cos 0. ~ Meaty. cos UG cos (wt+8_) JA} 


1 
— ed = 4. 
-) Wy x Vai + ve sin 0. + (Mote, sin wt i sin (Wt 6.) JA} 


1 


T= Oo 


e+ x Ya "easin 0. + Ya, "sa°°® 6. ae ee of, Ly, sin witty, cos wt J 


fd f 


- Vj AQ [¥g, sin(ut +6.) + 4 cos(ut +6,)]) 


iB | 
= Cree Xo + x 
A = e Doak nee 
c ac (r +r) 
Ss a 
= M + 
va ede tomtepe aa ©°° °s “deo 
Wy = -~ V._sin 6 _ F 
q1 Sa Ss qc 
x + x, a x" + xX. x +X. ay x" + xX. 
ose cl de p28 J 8. 
dc} x +x x +x x + x x +x 
Ss Ss Ss Ss 
t t 
x' ae x" rn x = x'' x - x" Te x on: x" 
-_ dc ee a Se dc see 
nace xT + xX © " x'+ x x + X . - x + x, 
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t 
= x! a T" 
wee x + Xo ic aS 


The preceding equations give the time dependence of trapped fluxes and 

torques following faults; they will be used later in developing calcula- 
tions of rotor heating following faults. However, the peak torque occurs 
during the subtransient period, and is the critical value for predicting 
failure of the mechanical structure. Simplified expressions for the sub- 


transient torques are given on the following pages. 


Torques During Subtransient Period 


Three-Phase Short Circuit from Open Circuit at t = 1/2 
Vows 1 


a fd fo 
i eee cos wt 
v2 v2 
a a 
a Bee cos wt ene = ~ ae at Wt = 7 (AS53) 


Three-Phase Short Circuit from Load at wt = 1/2 + 6 


V 
T= - —{ [M 
x 





= meet —_ ld 
fd- fo (x-x )I, sin(6+6)]cos wt + [(x-x )T cos (S+8) sin wt} 
(A-34) 
Line-Line Short Circuit from Open Circuit at t = wt 
2 
V sin wt (M,-I. ) 
a oO df fo 
T= - — sr Ma ele, cos Wt + a sin 2wt 
Va = Mateo : we for wt =z 71/2 
tT =- vee /x" cos wt + Vserex. sin 2wt (A-35) 
2 
max a 4 x" 
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Torques during Subtransient Period, (cont'd) 


Line-Line Short Circuit from Load at Wt 





De cos (w- 6.) 
tT = etx (Va, Sin wt + van cos Wt } 


Ve sin(wt, - 6) 
ia Soe ee ee cos wt - van sin wt] (A-36) 


x 
= 2 2 
+ it ’ i ! 2 + — 
Pee + sun +x)! Ya s¥q, COS 2 wt + Ug, ~ Yq, ) sin 2ue] 


x + oe = x" 
a 7 x + X Mateo * x+ x. ee od 7 





x - x" 
sen Teeter x ee sin 6. 





Closing Out of Phase by Angle ee at t = 0 


1 
i x" +x, ES sin G5 + eee cos 0. + Megl eg l¥y sinut + Va cos wt] 


wale way sin(wt+0 _) + Dee cos (wt +9_) ] 


a 


Voltages matched Vad = Mealeo 


van a Nes vee > ) 
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2 


sa 
t = et x {sin 661 ~ cos wt) + sin wt(l - cos 6.) ] (A-37) 


Rotor Torques 
The total torque acting upon the stator is given by the following 


equations. 


a Varta - Varta a 


To determine the distribution of the reaction rotor torque between the 
shield and field winding, the shield torque may be written as follows: 


Poerenier tg enka (A-39) 


If we use the constraint upon inductances and the definitions of the pre- 
vious section, Eq. (A-39) can be resolved into the difference of total 


torque and a component depending only on field current. 





The field winding torque is given by the second term above: 


tT. ® _ -k’ta iw (A-40) 
f Me f"q) 
The alternating field winding torques can then be found by solving for 
the induced field winding ac currents since Vee has no ac component. 
For all short circuits from open circuit, van = 0; therefore, Song O. 
The largest component of field winding fault torque will be alter- 


nating, since unidirectional torques are associated with resistance 


losses which are extremely small for a superconducting winding. 





In Eq. (A-18), the alternating component of field current as a function 
of alternating direct-axis current was derived: 


i 
x 
Mat dic 


i el ee 

| f | ie 
ac do ff " 2 
Y1+ (T 5) 


(A-41) 


The alternating field winding torque is given by the following equation: 


Use of the Lumped-Parameter Model Results 

The lumped-parameter models have yielded expressions for the per- 
unit currents, total torques, and field winding torques following short 
circuits. However, for use in calculating the heating and normal forces 
to which the shield is subjected, it is more convenient to write the 
armature currents in forms which identify the traveling waves of current 
which move forward, backward, and are stationary with respect to the 


Stator. This can be done by assuming first harmonic spatial distribution 


of the current. 


4(68,t) = i (t) cos 6 + ig (t) sin 6 (A-43) 


When the fault currents are written in this manner, the currents can 
be resolved into component waves which move relative to the stator. 
Current components with cos (ot + 6) move backward, relative to the 
Stator, at rrequency Ww, and backward relative to the rotor at angular fre- 


quency 2w ; they will be designated i, and referred to as "negative 
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sequence’’. Components stationary with respect to the stator move back- 
ward relative to the rotor at angular frequency Ww; they will be desig- 
nated i. and referred to as a "zero sequence". 

Current components with dependence cos (w - 6) will exist before the 
short circuit if the machine is initially loaded, but will increase in 
value following the short circuit; these current components are station- 
ary with respect to the rotor, and create the magnetic field which is 
initially excluded from the electrothermal shield. The increase in these 
currents over the value before the short circuit will be designated Ai; 


this corresponds to increase in positive sequence currents. 
Tnree-Phase Short Circuit 


From load at ee = 71/2 + 6 


ae Yay a: 
i(6,t) = TF cos (wot - 6) - a sin (wt - 8) =n sin 6 





’ 


ao) 
Ala, iad (hau) 


V 
; / a 
mm 


Line-to-Line Short Circuit 


From open circuit at t=t 


O 
’ ¥ ; 
i(6,t) = sh eaieoye- ai sin 0 - sok cos (ut + 6) 
e.. 2x'' - a me 2x" paid 2x" 
a 3 me 
Mat at pie 4 FE Be TE + ae 
Vasin@get s=u6) - 2 
een eee eM 
O x 
V 
a 
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From load at t = C 


l ox +x. 
i(6,t) 2x" +x.) xt ee _ van cos 6 )cos (wt - 6) 
2x" + x. 
= Sas ae + ve sin 6 )sin(wt oa 6) 
Wy 
- =} sin 6 
a ee +V__ cos 6_)cos(wt + 6) 
2 (x +x.) x" "dy Sane ge ae 
Xs 
= vee - ee sin 6 sin (wt + 6) 
l 2x" + x. 
Als, = 2x" + x, ( =o ay = Veg cos 65 = I. sin(6 + o”) 
(A-46) 
sie ¥B0 
Oo x 
i 2 2 2 
| (1%, +i, ) 


1 s 
a ZO" +x,) | x" Wa a aa rae Seb 


163 7 oe +x)‘ x an ee sin 55) 
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A-2 The Quasistatic Magnetic Field Model 


The armature currents which result when a generator is subjected to 
various fault conditions have been derived in section A-l-b. They were re- 
solved into traveling waves of current moving relative to the electrother- 
mal shield and a positive sequence component which appears to remain 
stationary. In order to determine the electromechanical effects upon the 
shield, a two-dimensional magnetic field model of the armature-electro- 
thermal shield will be developed with traveling waves of armature current 
as the source of magnetic field. 

Two types of effects are of interest in this model. One is a quasi- 
Steady-state effect produced by a traveling wave of armature current which 
decays slowly relative to its angular frequency, producing induced currents 
in the moving conducting shield. Heating of the shield results, and normal 
forces due to the interaction of the induced current with the armature 
field tend to produce ring vibration of the shell. The second effect re- 
sults from the sudden increase in the positive sequence component of arma- 
ture current. Currents are induced in the shield to keep its total flux 
linked constant, and exclude the demagnetizing effect of the armature 
currents. The currents induced in the shield decay with a time constant 
which depends upon shield parameters, but a large magnetic pressure is 
exerted upon the shield during this transient. 

In section A-2-a, the steady-state magnetic fields model will be 
derived for E.T. shields both thin and thick compared to skin depth. 

The results of these models will be used to predict heating and attenua- 
tion of the magnetic fields inside the shield. Experimental results are 


presented in Chapter V to verify the analytical results. In section 
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A-2-b, the transient-induced currents are investigated. Expressions for 
shield time constants are developed, and experimental verification is 
included. The magnitudes of the transient normal forces are also 


derived. 
A-2-a. Thin Shell Steady-State Mapnetic Field Model 


Figure A-5 rep- 
resents the cross- 
sectional view of the 
distributed armature 
current between radii 
Fi and Sy with the 
conducting E.T. shield 
at radius R. The 
boundary condition at 


radius T represents 


either infinite per- 





meability of infinite 
conductivities. The Figure A-5 

armature current is in the y- direction and has a sinusoidal spatial dis- 
tribution of n cycles and moving witn angular frequency WwW . The conduct- 
ing E.T. shield is rotating with angular velocity We The shield is assumed 
thin compared to the skin depth of the material at the frequency of the 
induced currents. The problem is to find: (1) the magnetic field outside 


the shield, (2) the induced currents and power dissipated in the shield, 
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and (3) the magnetic field inside the E.T. shield as a function of 
shield parameters. 
The problem is solved by replacing the distributed armature current 
by a current sheet. The effect of the distributed current is obtained by 
integrating the current sheet solution. 


a. rie? 


Definitions for an annular region 








with outer radius a and inner 


radius 8 are as follows: 


B = HH (r=0) 


B 


NMDrH R 


= UH (rsp) 
Ky = - Hy (r=a) 
Ka = Hy (r=8) 
The a and @ notation are 
replaced by a-b, c-d, e-f 
for the three annular regions 


in Fig. A6. 


Figure A-6 


The solutions for Laplace's equation for each of the annular regions 


a-b, c-d, e-f are combined according to the following boundary conditions: 


Boundary Conditions 


r=T BY = 0 for u> © Your 
or K? » O foro *+© ro>wT 


(wherever a double sign appears in the solution (+), 


the top sign applies for the iron (u * ©) outer shield, 
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and the bottom sign for the image (0 + ™) shield. 


r=S Bees 
ay r 
K +K =K 
Cc n 
r=R Bo = Bo 
r r 
K +K,= Asorw - 2 pf 
e d m n r 
r = 0 Ke = ( 


The solution is as follows, with Wom i — 


s\e 


The tangential magnetic field outside E.T. shield (r = R*) is" 


- Yl + (AouRw )? ( ~ hy a 
URW, 4(tan af tan OURW  ) 
Hae * eet e Bese a 
1+ (T Ww) ae 
sn xX 


The tangential magnetic field inside shell (r = R ) is: 


-1 
Hwee K J a j tan T,,4, 4 


1) At wa) ie 
x sn 


H is the magnitude of the magnetic field intensity which would exist 


w/o n 


at the shell surface if the shell were removed and the stator currents 


held constant. 





ie S.on 

x7 ©) aes ) IK, current sheet as (A-48) 
H = = 7 n+2 
w/o n a i 52-8 yy? a pom} S, ext?) 3 

2 2-n “o i = 72n n+2 n distri- 


buted stator 


current, as in Fig. A-5. 


a 
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Current induced in shield 
i 


= 
AouRw j(1 + tan = 
Kk. = ee ee Tn Wx oe (A-49) 
VEEL) Wile 
x Sn 


The time-average power dissipated in the shell is calculated from 


2 
IK. | 


PS “Lee The attenuation factor is the ratio of magnetic field 


inside the E.T. shield to magnetic field intensity without a field, but 


with the same armature current. 


Time-Average Power Dissipation (Thin E.T. Shield) 


2 i ale 

R 1 n m sn _ 2 
2 

[lt@)] ‘sn 1+(@-w)t_}? 


ene ton w/o n 


dn oO (A-50) 


where = 1 - R, 2n 
R - ree z ead turn 


Attenuation (Thin E.T. Shield 


[H,, | 
eee (A-51) 


AEC lH | Ww : Le 
w/o n {1 + Ems o Toa! } 
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A-2-b The Thick Shield, Steady-State Model 


The electrothermal shield can be considered thick if the skin depth 
of the induced currents is small compared to the shell thickness. This 
definition applied only for one frequency and the meaning of a thick shield 
is less apparent for transients. For steady-state power dissipation, the 
problem can be solved by simply replacing the electrothermal shield in 
Fig. A-5 by a solid cylinder of conductivity 0 + © and radius R. The 
boundary conditions are changed such that now Be = 0 at r=R. The result 
is familiar, and shows that the induced surface current is twice the mag- 


nitude of H 


acs The total power dissipation is then the same as would be 
W 


dissipated by the surface current distributed uniformly over one skin 


depth. The result is: 


Power Dissipation (Thick E.T. Shield) 
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The attenuation factor for a thick shield is much harder to obtain. It 
involves solution of the diffusion equation in cylindrical coordinates, 
which results in Bessel functions. The result can be simplified to the 


following equation: 
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The complete solution for the thick E.T. shield is presented in the 


following section. 
Solution of the Diffusion Equation for an Annular Region 


In section A-2-a, the solution for the induced currents in a thin 
conducting shield produced by armature currents was developed by combin- 
ing solutions to Laplace's equation for an annular region. For a thick 
conducting shield, the appropriate equation is the diffusion equation 
which can be derived from the Maxwell's equation neglecting displacement 
current and assuming that current density and electric field intensity 
are related by a constant conductivity. The resulting form of the 
diffusion equation is 

= Vl = (2 + v°V)H 
The probiem is to solve the diffusion equation for an annular region as 


*, K°, Kt which are 


in Fig. A-7, and find a relationship among B., B 
defined below, to represent the radial and tangential components of the 


magnetic field intensity at the outer and inner radii of the annular 


region. 


> 


WH (r=R,) 


oo 


y= Ho 


Cd > 


UH (r=R,) 


w~_> 
Oo 
i 


~ Hy (r=R,) 


> 
4. 


K = Hy (rs, ) 











—— —_ -« - Ew oe oOo =) om ~ 
® .& pp aaa s 7 
= = ———/?> oo ELF - 





153 


In two-dimensional cylindrical coordinates, the diffusion equation can 


be written for the radial component 


2 
ee ot) yh oe -(2.+52 2), 
uo \ r or or r2_ 08 ot xr o0/ r 


where v = Lov = Low or. 
Assume H = H_(r)ed (erne) solve for Hr) and obtain Ho from 
VeuH = QO. 
oH oH : eens 
2 — a as = 
rape ter yr (n° + ja°r dH 0 
q? = yoO(W - nw ) 
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The solution can be written as follows: 
— n at 
came Ot Cys qc 5 TK, ( v9 ar) 


Since we are interested primarily in the attenuation of the first spatial 


harmonic, we will let n = l. 
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The boundary conditions are given in the definitions of Bae Bs ss eRe. 
above. After much manipulation, a relationship among the field compon- 


ents at the edges of the annular region is obtained. 
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a' = (AF — BE) 


g' = (CF - DE) = - j 
ia =  <AH — BG) 

$' = (CH - DB) = - j 

w = AD - CB 


A = 73 aR, Tom vj aR) te oe al GEL) 
B= VjqRK (Yj aR) 5 Fo =§ K (79a R,) 
eee oe EG ak) 6» «Ga j J CF ak.) 
_ v3 aR, K, (Vj aR,) eae a K (v9 aR) 


The above expressions can be simplified by making assumptions about the 


magnitude of the skin depth of currents in the conducting annulus. 
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To find the attenuation of an alternating magnetic field by a thick 


shield, i.e., one for which A > 6, we assume a uniform magnetic field far 


==. @o w@& 
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from the thick shield alternating in magnitude at frequency W. This, 
of course, is equivalent to two counter-rotating magnetic fields rotating 


at frequency w. The solution for the ratio of field inside to outside is 


given by 
B 
a. - - 2 Bo 
By Rae (jw - 5°) (ju- a')+ 1 (A-52°) 


Further assumptions lead to increased 


simplification of w, a', and 6°. 
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The attenuation magnitude is given by 
as 
Lon pn (A-53) 
ie 


Thin Shield VY < 6 


In the thin shield case, the solution must reduce to that obtained 


in the previous section. 
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Equation (A-54) is identical with Eq. (A-51) obtained from the thin shell 
model. Experimental verification of Eqs. (A-53) and (A-54) is given in 
Chapter V. 
A-2-c Transient Magnetic Fields Model 

Thin Shield 

The time constant associated with decay of a sinusoidally distri- 
buted current sheet in a cylindrical shell of surface conductivity 0. 


can be calculated as follows: 


K = KA cos n@ att= Q 


K K_ cos n§@ 
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We assume time dependence of E, H to be 


-at 
= and substitute into the curl equation. 
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The value of 0. in terms of the bulk conductivity of a thin shell of 


thickness A is om = Ag 


HoAR 


an ; >= a (A-55) 


If the conducting shell is placed inside a boundary of 1 > * of ao > 


at r = T, the result will be 


1 R, 2™ 
ee 5 ULAR 3 + &) | 


It should be emphasized that the n refers to the harmonic of the initial 
current distribution. The longest time constant is for n= 1. The time 
constants associated with higher harmonics should not be confused with 
time constants associated with diffusion of the current through the thick 
ness of the material when a magnetic field is suddenly established out- 
side the shield. The diffusion time constants are derived in the next 
section. The result is that the first diffusion time constant is smaller 
than the decay time constant by a factor of the thickness-to-radius ratio. 

Thick Shield 

The terms "thick" and "thin" as applied to an electrothermal shield 
have significance only in terms of some steady-state frequency of induced 
eddy currents. However, it is clear that, for a shield whose thickness 
is larger than a small fraction of its radius, magnetic diffusion will 


result in time lags of the establishment of a magnetic field inside the 
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shield after a field has been extablished outside. This transient problem 
has been investigated [8, 18] and the solution for the time constants has 
been shown to reduce to the solution of the transcendental equation of 


Bessel functions 


Y, GR) ; ¥ , (O,R,? ee 
J, (o,R,) J, (a8, 


whe re R, and Ry are the inner and outer radii of the shield, and O. 


are related to the time constants of the field inside by 
tT = WW, (A-58) 
n om 


The magnetic field outside is uniform a long distance from the shield, 
and is established in a time short compared to the shield time constants. 


The resulting magnetic field intensity inside is then 


£0 _i 
gb 
= = nh 
H, (t) in) el 
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To determine the solutions for the Th approximations for the Bessel 


functions from Ref. [19] are used. 
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Using the above approximations, Eq. (A-57) reduces to 





15 1 
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The maximum value of aA is related to the longest shield time constant 
and the shield thickness. If we assume O << 1, and solve for To: the 
assumption can then be checked for typical dimensions of electrothermal 
shields. This has been done, and the value of a,A is usually less than 
.34. Using the assumption a4 < 1, we can expand the sines and cosines 


and obtain 





which gives the result for the time constant 


Aou R, i 


LA = Se AL ony (A-60) 


For A << Ro» this reduces to the result of the thin shield assumption ob- 
tained in the previous section. 

The higher order time constants are obtained from Eq. (A-59) by 
observing that, for larger n, the magnitudes of O. become smaller; i.e., 
a, are larger and the second and third terms can be neglected in Eq. (A-59). 


This gives 


7 (n- 1) 
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oO i 
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2 
n (n-1) 272 = 


Therefore, the first diffusion time constant (n= 2) is smaller than the 
L/R shield time constant (n = 1) by a factor less than the thickness-to- 
radius ratio. Experimental verification of these results is presented 


in Chapter V. 





169 


A-2-d Transient Normal Forces on a Conducting Shield 


When a magnetic field is established external to a conducting shield 
in a time short compared to the eddy current decay time of the shield, nor- 
mal forces are exerted upon the surface. These forces result from the 
interaction of the external tangential magnetic field and the surface cur- 
rents which are induced to prevent 
instantaneous change in the flux 
linked by the shield. In Fig. A-8, 
the external field which is uni- 
form far from the conducting shield 0 See 
has been established in a time short 
compared to the shield time const- a7 
ant. The currents in the shield are 


a oe ge 


e 
& 


uniformly distributed across the 

thickness ot a time long compared to 

the diffusion time constant. That Figure A-8 

is, if the field were established at t = 0, the observation is made at 
e=t., where roy <t. < = AGUR. It is clear from the figure that an 
inward radial force is being exerted upon the surface, the spatial de- 


pendence of which is cos’ 0. 


During transients in a generator, an increase in the positive 
Sequence component of armature current will produce a demagnetizing field 
which will initially be excluded from the electrothermal shield. As the 
induced currents in the shield decay, the armature produced field will 
penetrate the shield. During this transient, radially-directed forces are 


exerted which are stationary with respect to the shield. The trapped dc 


= Aas = 
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flux in the armature produces a traveling wave of normal stress upon the 
shield. Unbalanced faults such as line-to-line result in negative sequence 
magnetic field, which also produces traveling waves of normal stress. 

The magnitude and distribution of the normal stress is derived by 
calculating the distribution of magnetic field inside and outside the shell 
during the transient. The magnitude of the normal stress for a three- 
phase short circuit is independent of loading condition. The line-to- 
line fault is dependent upon loading, but the maximum normal stress which 
occurs one-half cycle after the fault is calculated for the open-circuit 
fault case. 

Normal Stress 
Three~Phase Fault from Open Circuit 
The stresses are independent of time 
of application of the fault. For simplifi- 
O 


cation, let t_ = 1/2 where t is time of 


application of fault, and t' equals time of 
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As. increase in direct-axis, positive sequence current following 
fault 

= increase in direct-axis, positive sequence current for time t 
such that Ty << t' << Ty 

The normal stress given by Eq. (A-61) is made up of a component indepen- 

dent of angular position, and one having a cos 2y dependence. Since the 

shell is much less resistant to the cos 2y stress, its maximum value is 


important. This occurs for wit = , i.e., one-half cycle after the 


fault. The second spatial harmonic normal stress is given by 
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A-3 The Thin Shell Mechanical Model 

During transients following machine faults, large normal forces are 
exerted upon the electrothermal shield. These forces consist of a de 
or stationary force relative to the shield distributed with a cos” @ 
angular dependence, and a traveling wave or alternating normal force. The 
magnitudes of these forces were derived in the previous sections. In order 
to predict the mechanical motion of the shell when subjected to these nor- 
mal forces, the thin shell elasticity equations are solved for the given 
loading. Steady-state deflections are calculated as a function of load- 
ing, and ring natural frequencies are derived by the energy method. 


The Thin Shell with Hinged Ends Radially Supported 





The shell is loaded as shown in Fig. 
Ac? . 
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m=0 n=0 
The loading due to the armature-produced 
magnetic flux occurs primarily for n = l 


and n = Q,1 which can be written as 


2 TX 
COC = OO cos sin -—— 
r oO o L 


The distribution of loading was shown in 


section A-2-c to have the cos*  depen- 





dence. The loading distribution along the 
axial length'of the shield depends upon Fig. A-9 
the armature active length add end turns 


relative to the end supports of the E.T. shield. Fourier analysis of the 
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flux distributions for armatures of the machines being built, and the 
geometries of larger machines, has shown that the sin 7x/2 term very 
closely represents the axial flux distribution. Higher harmonics in the 
axial distribution are very small relative to the fundamental. 


The partial differential equations which govern the deflection of a 


thin shell, when loaded by a normal stress on the surface, are given in 





Ref. [15] and reproduced here: 

fee OW tu ley 
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wnere u, v, W are the axial, azimuthal, and radial deflections, res- 
pectively, and G. is the radial load which is a function of x and ®@. 
The primary loading due to magnetic pressure has a cos” @ dependence, as 


shown in the previous section. The loading is then assumed 


1 1 ® 
Te 7 COs $)sin 


TX 


q (A-65) 
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Assume solutions of the form 


u = cos | > A, cos no axial deflection 
v = sin - z B sin nd azimuthal deflection , 
w = sin radial deflection 


a Cc, cos nd 
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The shell is assumed to be hinged, with radial support at the ends x = 0, 


but is free to move axially. 


After substitution of the assumed solutions into the equations, ap- 
plication of approximations concerning relative dimensions, and solution 


for the coefficients, the resulting radial deflection is 


TT R* Do 4 R eee 
ar oT ol 2 oe (+) ae == cos 2 ¢ (A-66) 
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2 R = shield radius 

oe > Tl 

+ & = length between shield and support 
4 = shield thickness 
E = Young's modulus 


Vv = Poisson's ratio 


The ratio of Cy to a gives the relative magnitude of the uniform radial 


deflection to the second harmonic deflection. 





T — voor ie 
FIG) + 3) 
WL - vp 


2 
o 


2 
which is clearly less than 1. Even for the limiting case of 2 = TR, 


the uniform deflection is only about 1/16 the second harmonic amplitude. 


Shell Mechanical Natural Frequencies 


The natural frequency of ring vibrations having radial deflections 


of the lowest mode 


TX 
Wir x = w cos Wt cos 2¢ sin — 
(r5,x) ee ¢ sin | 


can be determined by equating the maximum kinetic energy of the motion 


to the maximum energy of bending. 





166 


To calculate the potential energy stored in a ring mode deflection, 

1 1 

~P6 +=P6 + 
2) 1 A eee 


Paes where P is the load in the 1 direction and oe is the deflection in 


we use Castigliano's theorem [16] which states that U = 


the 1 direction, etc. 
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6. Ww) = dr cos 26 sin 7 
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The kinetic energy maximum is calculated by assuming deflection time 


dependence as 





6 = Ww cos 26 sin — sin wt 
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wa Tl 
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The contribution of the axial deflections to kinetic energy is neglected, 


since for a long shield, they are much smaller than radial deflections: 


w~ 

t= 

i! 
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The energy change due to motion in centrifugal force field of a shield 


rotating at frequency O ([17} is: 
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ae eae o A R* w we 
Cc 2 m m max 


If the sum of the maximum potential and kinetic energies are equated, the 
natural frequency of the shield is obtained for the assumed mode of vi- 


bration: 


w* = 2.4w*+ w ? 
c m 1 
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WO. * natural frequency of rotating shield 
_ angular frequency of rotation 


pe density of shield material (kg/m?) . 
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